AD-A256  681 

*  . .  .....  nm  attl  nil  ail 


PL-TR-92-2149 


% 


THE  RECEIVER  STRUCTURE  BENEATH  THE 
CHINESE  DIGITAL  SEISMOGRAPH  NETWORK 
(CDSN)  STATIONS:  PRELIMINARY  RESULTS 


S .  Mangino 
J.  Ebel 


!  Weston  Observatory1^ 

"-Department  of  Geology  and  Geophysics 
Boston  College 
Weston,  MA  02193 


30  April  1992 


Scientific  Report  Mo.  l 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


PHILLIPS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 

HANSCOM  AIR  FORCE  BASE,  MASSACHUSETTS  01731-5000 


92-23675 


92  8  25  108 


1 


The  views  and  conclusions  contained  in  this  document  are  those  of  the  authors 
and  should  not  be  interpreted  as  representing  the  official  policies,  either 
expressed  or  implied,  of  the  Air  Force  or  the  U.S.  Government. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


ES  F.  LEWKOWICZ 
Contract  Manager 
>olid  Earth  Geophysics  Branch 
Earth  Sciences  Division 


Solid  Earth  Geophysics  Branch 
Earth  Sciencs  Division 


D^ML^^E^K^DT^rector 
Earth  Sciences  Division 


This  document  has  been  reviewed  by  the  ESD  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 

Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center.  All  others  should  apply  to  the  National  Technical 
Information  Service. 

If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization,  please  notify 
PL/IMA,  Hanscom  AFB  MA  01731-5000.  This  will  assist  us  in  maintaining  a 
current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on 
a  specific  document  requires  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No  0704-0188 

Public  reoortmg  burden  for  thi*  collection  of  .^formation  >s  estimate  to  werag*  *  hour  per  re* ponse.  including  time  »c r  reviewing  instructions,  snart-  -g  e*  sting  data  sources, 
gathering  and  maintaining  the  data  needed,  and  comoietmg  ang  reviewing  the  collection  o*  information  *>«*n d  comments  regarding  this  burden  est-^ate  G *  an^  other  aspect  of  this 
icliection  of  information,  including  suggestions  for  reducing  this  Durden  to  Washington  Heaaauarters  S erv  res.  Directorate  for  information  0:vm-.i: ns  ■*  a»oc;rts.  t2 IS  iefferson 
Davis  Highway,  Suite  1204.  Arlington.  VA  22232-4  302  and  to  tn*  Offee  of  M  i«*g*m*nt  and  Budget.  Paperwork  Reduction  Project  (07C4-0 ' 83).  Washing:.-  r  Z-C  20SC3 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  ANI 

30  ADril  1992  Scientific  Re 

3  OATES  COVERED 

port  No.  1 

4.  TITLE  AND  SUBTITLE 

The  Receiver  Structure  Beneath  the  Chinese  Digital 
Seismograph  Network  (CDSN)  Stations:  Preliminary  Results 

5.  FUNDING  NUMBERS 

PE  61102F 

PR  2309  TAG2  WU  AW 
Contract  F19628-9 1 -K0009 

6.  AUTHOR(S) 

S.  Mangino  and  J.  Ebel 

5 

j  7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

!  Weston  Observatory 

i  Department  of  Geology  and  Geophysics 
!  Boston  College 

Weston,  MA  02193 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING /MONITORING  AGENCY  NAM£(S)  AND  AOORESS(ES) 

Phillips  Laboratory 

|  Hanscom  AFB ,  MA  01731-5000 

1 

\ 

|  Contract  Manager:  James  Lewkowicz/GPEH 

10.  SPONSORING  MONITORING 

AGENCY  REPORT  NUMBER 

PL-TR-92-2149 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release; 

Distribution  unlimited 

12b.  distribution  code 

1  13.  ABSTRACT  (Maximum  200  words)  ] 

|  This  report  provides  preliminary  results  of  an  investigation  of  the  crust  and  upper-mantle  seismic 
velocity  structure  beneath  the  China  Digital  Seismograph  Network  (CDSN).  Receiver  functions 
j  obtained  from  over  100  teleseismic  events  were  stacked  by  common  backazimuth  and  epicentral 
distance  and  inverted  for  the  velocity  structure  beneath  5  CDSN  stations.  The  inversion  results  show 
that  stations  HIA  and  BJI  have  well  constrained  crustal  models  of  about  40-km  thickness  with 
gradational  Mohos.  At  WMQ  there  is  evidence  for  shallow  dipping  structure  and  the  inversion 
results  indicate  a  smooth  positive  gradient  in  the  crust  with  a  transitional  crust-mantle  boundary. 
Model  estimates  for  KMI  and  LZH  are  poorly  constrained  by  the  data  analyzed  thus  far.  Synthetic 
seismograms  of  5  CDSN  crustal  models  for  the  Soviet  JVE  nuclear  test  indicate  that  the  receiver 
structure  most  strongly  affects  the  radial  component  of  the  recorded  P  wave. 


14.  SU81ECT  TERMS 

Seismology  Crustal  Structure  Upper  Mantle  Structure 

China  Receiver  Function  Regional  phases  Teleseismic 

IS.  NUMBER  OF  PAGES 

108 

16.  PRICE  COOE 

18  SECURITY  CLASSIFICATION 

19  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

_ Unclassified _ 

Unc  lass i f i ed 

Unc  lassif ied 

SAR 

NSN  7S40  0'  280  SSOO  Standard  -orm  298  (Rev  2  89) 


Table  of  Contents 


List  of  Figures  .  iv 

List  of  Tables  .  v 

1.0  Introduction  .  1 

2.0  Geologic  of  China  and  crustal  structure  at  six  CDSN  stations .  4 


2.1 

BJI 

2.2 

KMI 

2.3 

LZH 

2.4 

MIA 

2.5 

MDJ 

2.6 

WMQ 

3.0  Receiver  Functions  .  18 

4.0  The  Data .  23 

5.0  CDSN  Receiver  Function  Observations . . .  27 


5.1 

BJI 

5.2 

KMI 

5.3 

LZH 

5.4 

HIA 

5.5 

MDJ 

5.6 

WMQ 

6.0  Receiver  Function  Inversion  .  41 

7.0  Preliminary  Results .  44 


7.1  BJI 

7.2  KMI 

7.3  LZH 

7.4  11IA 

7.5  MDJ 

7.6  WMQ 


8.0  Discussion  .  51 

8.1  Implications  for  monitoring 

9.0  References  .  56 


10.0 


Appendix 


Accesion  For  ^ 

NTIS 

CRA&I 

a 

0TFC 

TAB 

□ 

Unannounced 

Li 

Justification 

By . 

Dist.  ibution  / 

i  i  i 

— 

Availability  Cnc  es 

Avail  <. 

»!,  Ci 

Dist 

Sp« 

C.ul 

DTIC  QUALTT  If  Ill  OBJECTED  1 

List  of  Figures 


Figure  la  Physiographic  map  of  China  and  surrounding  region  with  the  most  prominent 

geologic  features .  10 

Figure  lb  The  provinces  of  China,  from  Hsu  et  al.  (1990) .  1 1 

Figure  2  Contour  map  of  Moho  depth  across  mainland  China  from  Yuan  et  al.  (1986).  ...  12 

Figure  3  Estimated  Pn  velocity  distribution  across  mainland  China  from  Xingxin  (1988)  13 

Figure  4  Location  of  seismic  refraction  profiles  east  of  station  BJI  and  the  corresponding 

modeling  estimates  from  Changquan  and  Shixu  (1986) .  14 

Figure  5  Estimated  P  wave  velocity  structure  near  KMI  from  Qizhong  et  al.  (1985) .  15 

Figure  6  Contour  map  of  Moho  depth  for  the  Lanzhou  region  (top)  from  Xingxin  (1988) 

and  local  distribution  of  active  faults  (bottom)  from  Burchfiel  et  al.  (1991 ) .  16 

Fiuure  7  Crustal  structure  inferred  by  surface  wave  modeling  across  northern  China  from 

Rui  et  al.  (1981) .  17 

Figure  8  Characteristic  teleseismic  event  from  the  Fiji  area  recorded  at  station  BJI .  21 

Figure  9  Receiver  function  ray  diagram  and  corresponding  radial  receiver  function  (top). 

Synthetic  response  for  30  to  90  degrees  epicentral  distance  (below) .  22 

Figure  10  Distribution  of  CDS  N  data  for  calender  years  1986-199! .  25 

Figure  1 1  Flow  chart  of  procedure  used  in  this  report .  26 

Figure  12  Receiver  functions  of  station  BJI:  BJI-45-NE,  BJI-74-SE  and  BJI-56-NW .  32 

Figure  13  Receiver  functions  of  station  KMI:  KMI-63-NE,  KMI-56-SE,  KMI-64-SW  and 

KMI-46-NW .  33 

Figure  14  Receiver  functions  of  station  LZH:  LZH-55-NE,  LZH-46-SE  and  LZH-43-W.  34 

Figure  15  LZ11-SE  radial  vs  tangential  stacks.  The  apparent  moveout  of  a  PpPhs  multiple 
is  consistent  with  the  response  to  a  common  depth  boundary  with  increasing 
epicentral  distance .  35 

Fit: ure  16  Receiver  functions  of  station  HIA:  HIA-50-NE,  11IA-62-SE,  HIA-34-SW  and 

H1A-50-NW .  36 

Figure  17  HIA-SE  radial  vs  tangential  stacks  with  an  azimuthal  separation  of  45°.  In 
comparison  to  LZH  the  tangential  energy  is  minimal  at  HIA.  Arrivals  at  5-5. 3s 
an  at  12-  14s  are  consistent  with  a  Ps  and  PpPh-s  multiple  converted  phases  from 
the  crust  mantle  boundary  as  the  ray  parameter  decreases .  37 

Figure  18  Receiver  functions  of  station  WMQ:  WMQ-63-NE,  WMQ-63-SE,  WMQ-50- 
SW  and  WMQ-38-NW . 

1  V 


38 


Figure  19a  Variation  of  peak  amplitude  first  arrival  time  as  a  function  of  backazimuth  for 
the  most  stable  WMQ  receiver  functions.  The  pattern  suggests  a  dipping 
structure  (Owens  and  Crosson,  1988b) .  39 

Figure  19b  Tangential  receiver  functions  as  a  function  of  backazimuth  for  WMQ.  Polarity 
reversals  (denoted  by  arrows)  occur  along  a  line  parallel  to  the  dip  direction 
(Langston,  1977) .  40 

Figure  20  Suite  of  48  starting  models  generated  for  station  BJI  .  43 

Figure  21  Preliminarily  modeling  results  for  station  BJI .  46 

Figure  22  Preliminarily  model  synthetics  for  station  KMI .  47 

Figure  23  Preliminarily  model  synthetics  for  station  LZH .  48 

Figure  24  Preliminarily  modeling  results  for  station  HIA .  49 

Figure  25  Preliminarily  modeling  results  for  station  WMQ .  50 

Figure  26  Great  circle  arc  paths  from  Semipalatinsk  to  CDSN  stations .  53 

Figure  27a  The  first  20  seconds  of  the  Soviet  JVE2  nuclear  test  as  recorded  by  the  CDSN..  54 

Figure  27b  Synthetic  seismograms  convolved  with  a  source  time  function  appropriate  for 

the  J  VF2  and  instrument  response  for  each  CDSN  station .  55 

List  of  Tables 

Table  1  Event  listing  and  distribution  for  station  BJI .  59 

Table  2  Event  listing  and  distribution  for  station  KMI .  60 

Table  3  Event  listing  and  distribution  for  station  LZH .  61 

Table  4  Event  listing  and  distribution  for  station  HIA .  63 

Table  5  Event  listing  and  distribution  for  station  MDJ .  64 

Table  6  Event  listing  and  distribution  for  station  WMQ .  65 

Table  7  CDSN  station  locations .  67 

Table  8  CDSN  sensitivity  STS-l-BB  channels .  68 


V 


The  Receiver  Structure  Beneath  the  China  Digital  Seismograph  Network 
(CDSN)  Stations:  Preliminary  Results 


Stephen  Manf>ino  and  John  Ehei 


Weston  Observatory 
Boston  College 
Weston,  MA  02193 


Summary 

This  study  provides  preliminary  model  estimates  of  the  crust  and  upper  mantle  velocity 
structure  in  China  beneath  the  Chinese  Digital  Seismograph  Network.  The  one-dimensional 
models  are  derived  from  receiver  functions  obtained  from  teleseismic  data  recorded  at  each 
station.  Individual  models  are  used  to  investigate  the  effects  of  the  receiver  structure  on  regional 
waveforms  recorded  at  each  station. 

Key  Words:  receiver  function,  teleseismic,  inversion,  shear  velocity,  crust,  upper  mantle.  China. 

l.o  Introduction 

A  renewed  interest  in  furthering  restrictions  on  testing  and  proliferation  of  nuclear  weapons 
has  motivated  the  seismological  community  to  develop  effective  regional  and  teleseismic 
methods  of  detection  and  discrimination  between  small  earthquakes  and  low  threshold  nuclear 
explosions.  In  particular,  the  availability  of  high  quality  data  from  seismograph  stations  installed 
in  China  and  the  Russian  States  has  opened  up  research  efforts  to  improve  current  levels  of 
effective  detection  and  discrimination  of  events  recorded  on  the  Eurasian  continent. 

Most  current  regional  explosion  yield  estimation  and  discrimination  methods  are  based  on 
spectral  analysis  of  seismic  data,  however  the  spectra  of  particular  regional  seismic  phases  such 
as  Pg,  Pn,  Sn,  and  Eg  are  affected  by  attenuation  and  scattering  near  the  source,  along  the  crust- 
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mantle  path,  and  near  the  receiver.  As  a  consequence,  for  a  given  region  correction  terms  are 
necessary  prior  to  evaluating  data  from  a  given  station.  One  way  to  estimate  the  crustal  velocity 
structure  is  by  active,  dense  seismic  refraction  and  reflection  surveys,  but  unfortunately  such 
costly  efforts  are  often  not  practical  or  feasible.  An  alternative  method  is  through  the  detailed 
analysis  of  passively  recorded  teleseismic  data.  The  purpose  of  this  work  is  to  determine  the 
crustal  response  at  several  stations  in  China  using  teleseismic  data  in  order  to  better  understand 
earthquake  and  explosion  waveforms  recorded  by  these  stations. 

Teleseismic  body  waveforms  provide  a  wealth  of  information  on  the  velocity  structure 
beneath  three-component  seismic  stations  (for  example,  Archambeau  et  al.  1966;  Basham  and 
Ellis  1969;  Kurita  1969;  Bakun  1971).  Time  domain  modeling  of  broadband  teleseismic  P- 
waveforms  to  determine  the  crust  and  upper  mantle  velocity  structure  beneath  three  component 
seismic  stations  is  well  established  in  seismology  (Burdick  and  Langston  1977;  Langston  1977b; 
Langston  1979;  Langston  1981;  Owens  et  al.  1984,  1987,  1988a,  1988b;  Ammon  et  al.  1990; 
Mangino  et  al.  1992).  In  this  report  we  describe  the  response  and  present  preliminary  models  of 
the  seismic  velocity  structure  beneath  five  China  Digital  Seismograph  Network  (CDSN)  stations. 
Operated  by  the  State  Seismological  Bureau  Peoples  Republic  of  China  in  cooperation  with  the 
Albuquerque  Seismological  Laboratory  US  Geological  Survey,  the  CDSN  stations  record  high 
quality  broadband  digital  data  in  continuous  and  trigger  mode.  This  network  currently  represents 
12 %  (as  of  March  26,  1992)  of  all  stations  consistently  providing  data  to  the  Global  Digital 
Seismic  Network  (GDSN).  Across  mainland  China  from  east  to  west  shown  in  Figure  la  the 
CDSN  stations  are:  Mudanjiang,  Heilongjiang  province  (MDJ);  Hailar,  Neimenggu  Province 
( H I A ) ;  Beijing,  Baijatuan  (BJI):  Kunming,  Yunnan  Province  (KM1);  Lanzhou,  Gansu  Province 
(LZIli;  and  Urumqi,  Xinjiang  Province  (WMQ). 

Sections  1  and  2  summarize  current  interpretations  of  the  tectonic  history  and  describe  some 
a  priori  crustal  information  which  was  gathered  in  the  vicinity  of  each  station.  We  use  this 
information  and  forward  modeling  to  construct  initial  models  for  the  time  domain  inversion 
procedure  (Ammon  et  al.  1990)  used  in  this  report.  Section  3  summarizes  the  receiver  function 


modeling  technique  and  in  Section  4  we  describe  the  teleseismic  data  obtained  from  the  archives. 
Section  5  describes  the  receiver  functions  at  each  CDSN  station.  The  responses  at  each  station 
ranges  from  simple  to  complex,  and  the  degree  of  complexity  seems  to  correlate  with  the  amount 
of  geologic  deformation  experienced  at  each  site.  Section  6  presents  a  summary  of  the 
teleseismic  P-waveform  time  domain  inversion  procedure  given  in  Ammon  et  al.  (1990).  Section 
7  presents  preliminary  estimates  of  the  crustal  velocity  structure  beneath  each  CDSN  station.  To 
date  our  focus  has  been  on  matching  only  the  most  prominent  arrivals  in  the  data  within  ±1 
standard  deviation  bounds  obtained  from  the  variance  of  the  stacked  data.  Section  8  assesses  the 
implications  for  monitoring  earthquake  and  explosion  sources  at  these  stations. 
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2.0  Geologic  Overview  of  China 

The  following  sections  are  intended  to  summarize  the  current  view  of  the  tectonic  history  of 
mainland  China,  in  order  to  build  a  framework  for  interpreting  the  observed  receiver  functions 
and  inferred  crustal  models  beneath  each  CDSN  station.  One  problem  encountered  in  preparing 
this  summary  was  the  lack  of  consistent  nomenclature,  ubiquitous  among  English  language  and 
translations  of  Chinese  papers,  in  reference  to  the  same  crustal  features  and  physiographic 
province.  From  previous  work  we  can  obtain  estimates  of  the  average  crustal  P-wave  velocity. 
However,  estimates  of  total  crustal  thickness  vary  by  as  much  as  8  km  and  upper  mantle  P 
velocity  vary  by  as  much  as  0.5  km/s  from  paper  to  paper  for  some  stations.  Overall,  Acta 
Seismologica  Sinica  edited  by  the  Seismological  Society  of  China  was  a  most  useful  source  of 
information.  For  additional  information  on  the  tectonic  history  and  current  topics  of  research 
refer  to:  Yuan  et  al.  (1986);  Hsu  et  al.  (1990);  Zhaoet  al.  ( 1990);  Gilder  et  al.  (1991);  Shangyou 
( 1991)  and  Birchfiel  et  al.  (1991). 

China  may  be  considered  a  composite  continent  of  several  paleo-plates  delineated  by  major 

suture  zones  and  fold  belts.  The  “nucleus”  of  present  day  China  centers  on  the  Tarim  and  Sino- 

Korean  (North  China  Block)  cratons  onto  which  several  additional  land  masses  have  accreted 

(Gilder  et  al.  1991).  Figure  la  shows  the  principal  tectonic  subdivisions,  major  suture  zones  and 

prominent  crustal  features  of  mainland  China.  Shown  in  Figure  lb  are  the  modern  political 

boundaries.  Accretion  to  the  northwest  of  the  nucleus  are  marked  by  the  Tien  Shan-Tumen  and 

Junggar-Xingan  fold  belt  and  the  Late  Paleozoic  Junggar-1  legen  suture  to  the  Siberian  craton 

(Gilder  et  al.  199 1 ).  Southeast  and  southwest  of  the  Tarim  and  North  China  Block  several  stages 

of  accretion  have  been  identified,  but  the  relative  timing  and  emplacement  of  individual 

microplates  are  currently  a  focus  of  debate  (Hsu  et  al.  1990;  Huang  and  Opdyke  1991;  Gilder  et 

al.  1991:  and  others).  To  the  southeast  accretionary  events  are  represented  by  the  east-west 

trending  Kunlun-Quinlin  fold  belt  and  by  the  Yangtze  Hunan  Deformed  Belt  which  is  also 

known  as  the  South  China  Block  (Gilder et  al.  1991).  Southwest  of  the  Kunlun-Quinlin  fold  belt 

are  the  Tangla  Deformed  Belt  and  Lhasa  block,  and  since  the  Mesozioc  the  collision  between  the 
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Eurasian  and  Indian  Plate  has  formed  the  Himalayan  Mountains  and  Tibetan  Plateau  (Gilder  et 
al.  1991;  Shangyou  1991). 

Between  100-107°E  a  roughly  north-south  seismogenic  zone  divides  the  whole  country  into 
two  parts.  In  general,  the  eastern  half  of  China  is  affected  mainly  by  circum-Pacific  tectonics 
and  associated  back-arc  extension,  while  the  western  part  is  primarily  affected  by  Tethys- 
Himalaya  compressional  tectonics  (Yuan  et  al.  1986).  China  has  an  extensive  historical  data 
base  of  earthquakes  due  to  the  early  interest  in  seismological  study  in  132  AD  by  Chang  Heng 
and  the  invention  of  the  seismoscope,  the  world’s  first  seismometer.  Fu  ( 1991)  tabulates  Chinese 
earthquake  data  between  19(H)- 1980  and  statistically  shows  that  the  seismicity  of  eastern  and 
western  China  is  dominated  by  shallow  focus  intraplate  earthquakes.  These  earthquakes  occur  7 
times  more  frequently  in  western  China  compared  to  eastern  China  and  have  a  corresponding 
energy  release  ratio  of  23: 1  over  the  same  time  period. 

Broad-scale  information  on  crustal  thickness  and  upper  mantle  velocity  across  China  is 
shown  in  Figures  2-3  from  geophysical  maps  given  in  Yuan  et  al.  (1986)  and  Xingxin  (1988).  In 
Figure  2  the  estimated  crustal  thickness  along  the  eastern  coast  of  China  is  34-36  km  and 
increases  to  44-46  km  in  central  China  near  Chengdu.  West  of  Chengdu  the  Moho  depth  is 
inferred  to  be  as  much  as  70  km  beneath  the  Tibetan  Plateau.  These  results  are  consistent  with 
those  of  Zhao  et  al.  (1991)  who  infers  crustal  thickness  from  modeling  both  regional  Love  and 
Pni  waves.  Wu  (1990)  conducted  an  extensive  analysis  of  local,  regional  and  teleseismic  data 
recorded  at  Chinese  Kirnos  and  CDSN  stations.  For  eastern  China,  Wu  indicates  Pn  and  Sn 
arrivals  from  events  within  10°  are  generally  weak  and  unclear,  but  in  western  China  in  the 
Xinjiang  and  Tibet  regions,  Pn  and  Sn  may  be  quite  prominent.  In  reference  to  Figure  3, 
relatively  Large  Pn  velocities  are  inferred  for  the  Junggar  and  Tarim  Basins  and  for  the  western 
third  of  the  Sino- Korean  Block.  Much  of  this  map  is  speculative,  so  that  caution  must  be  taken 
in  extracting  absolute  values.  However,  the  general  trends  are  informative  and  we  make  use  of 
those  trends  in  this  study. 
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2.1  BKIJING  (BJI) 


Station  BJI  is  located  in  the  Sino-Korean  block  in  the  northern  end  of  the  Jizhong  Depression 
approximately  20  km  north  of  Beijing.  The  station  is  installed  in  a  hilly,  mostly  agricultural 
region  in  a  subsurface  vault  at  3  meters  depth  on  piers  resting,  on  gravel  (Peterson  and  others — 
Kexin  and  others,  19X7).  Topography  in  the  vicinity  of  the  station  is  minimal.  The  Jizhong 
Basin  is  characterized  by  several  NNE-SSW  trending  linear  basement  high  and  low  graben 
structures  composed  of  pre-Cenozoic  rocks  overlain  by  recent  fluvial  and  lacustrine  sediments 
(Zhao  and  Windley,  1990).  Gilder  et  al.  (1991)  discuss  two  stages  of  rifting  for  the  Bohia  Gulf 
Region,  the  most  recent  being  Late  Cretaceous-Eocene  until  present  with  a  heat  How  average  of 
1.6  IIEII  (1  HFU  -1  peal  env-s'1).  Since  the  1976  Tangshan  earthquake,  one  the  most 
devastating  earthquakes  (Mb=  7.8 )  in  world  history,  several  studies  on  the  crustal  structure  in  the 
vicinity  of  BJI  have  been  published.  Changquan  and  Shixu  (1986)  present  and  model  seismic 
refraction  data  from  two  refraction  profiles  almost  200  km  east  of  BJI.  In  reference  to  Figure  4, 
crustal  P-wavc  velocities  range  from  6. 1-6.4  knt/s  with  a  prominent  low  velocity  zone  located  in 
the  mid-  to  lower  crust,  and  crustal  thickness  increases  from  SE  towaru  NW  over  a  transitional 
crust  mantle  boundary.  Exclusive  of  the  low  velocity  zone  I  Iellinger  et  al.  (1985)  infer  a  crustal 
thickness  of  35-38  km  southeast  of  BJI.  Shedlock  and  Roecker  (1987)  infer  an  upper-crust  P- 
uave  velocity  at  6.02  km/s  between  0-15  km,  6.5  km/s  between  15-35  km  and  an  upper  mantle 
velocity  of  7.9  km/s  east  of  BJI.  They  also  infer  a  lower  crustal  low  velocity  zone  between 
Beijing  and  Tangshan  at  20-35  km  depth  as  well  as  Vp/Vs  ratios  of  1.73-1.77  in  the  crust  and  1.8 
in  the  upper  mantle.  Yuan  et  al  ( 1986)  estimate  a  36  km  thick  crust  beneath  BJI. 

2.2  Kl  wiim;  (KM I) 

Station  KMI  is  located  along  the  southwestern  edge  of  the  Yangtze  or  South  China  Block  in 
the  southern  Yunnan  province.  The  seismometers  are  installed  in  a  30  meter  long  limestone 
tunnel  27  meters  deep  and  are  within  5  km  of  some  major  industrial  plants  (Peterson  and 
others  Kexin  ami  others,  1987).  Topography  around  KMI  is  moderate  with  elevations  between 

f. 


1500-3000  meters.  Paleomagnetic  evidence  given  in  Fang  et  al  (19X9)  suggests  that  the  Yunnan 
province  is  composed  of  at  least  two  accretionary  blocks  (eastern  and  western)  and  that  either  the 
Red  (Yuan)  River  or  Lancang  (Mekong)  River  fault  zones,  -200  and  -300  km  west  of  KMI 
respectively,  represent  a  major  paleo-suture  of  the  Yunnan  and  South  China  Block.  Baiji  and 
Ning  (19X7)  delineate  a  ‘central’  Yunnan  block  bounded  by  the  Red  River  fault  to  the  south,  the 
Xiaojiang  fault  to  the  east,  and  the  Jianchuan-Lijiang  fault  to  the  north  shown  in  Figure  5.  Using 
a  network  of  14  seismic  stations  and  KMI  as  the  ‘standard  station’  Baiji  and  Ning  (19X7)  use 
teleseismic  P-wave  residuals  to  infer  a  higher  velocity  upper  mantle  beneath  the  central  block 
(0.45  km/s  faster)  than  the  surrounding  region.  The  result .  of  a  deep  seismic  sounding  profile 
recorded  along  the  main  road  from  Malong,  about  100  km  NE  of  KMI,  through  Kunming  to 
Simao  -300  km  southwest  of  KMI  is  discussed  in  Qizhong  et  al.  (19X5).  Their  model  estimate 
east  of  KMI  (Figure  5)  indicates  a  45  km  thick  5  layer  crust  with  an  average  velocity  of  6.3  km/s 
over  an  upper  mantle  of  -X.2  km/s.  Yuan  et  al.  (19X6)  estimate  a  total  crustal  thickness  beneath 
KMI  between  4X-50  km. 

2.3  Lanzhou  (LZH) 

Station  LZH  is  located  on  the  northeastern  margin  of  the  Tibetan  Plateau  along  the  Qinling 

fold  belt  in  the  Gansu  Province.  The  seismometers  are  located  on  the  outskirts  of  Lanzhou  City 

in  a  buried  surface  vault  built  on  loess  deposits  (Peterson  and  others — Kexin  and  others,  1987). 

Topography  in  the  vicinity  of  LZH  is  characterized  by  smooth  rounded  slopes  with  narrow,  steep 

walled  modern  stream  valleys  at  elevations  of  2(XX)-3()(X)  meters  (Burchfiel  et  al.  1991).  Several 

large  (Mb  >  6)  earthquakes  have  occurred  in  this  area  with  the  largest  being  the  1920  Haiyuan 

event  that  had  an  associated  220  km  surface  rupture  (Burchfiel  et  al.  1991).  Shown  in  Figure  6 

(top)  are  the  Quaternary  active  faults  and  Cenozoic  basins  about  Lanzhou.  Motion  along  these 

faults  suggest  the  Tibetan  Plateau  is  extending  toward  the  northeast.  Lanzhou  overlies  a  zone  of 

high  electrical  conductivity  in  the  upper  mantle  and  local  hot  springs  provide  source  regions  for 

high-temperature  wells  about  Lanzhou  (Xingxin  19XX).  Shown  by  contour  plot  in  Figure  6  from 
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Xingxin  (1988)  station  LZH  appears  to  overly  a  southwest  dipping  Moho  at  54  km  depth.  Liang 
and  Li  (1991)  infer  51-52  km  thick  crust  using  gravity  data  and  show  a  similar  dipping  nature  of 
the  crust-mantle  boundary  beneath  LZH.  Yuan  et  al.  (1986)  infers  a  50-54  km  thick  crust  (Figure 

2)  beneath  LZI I. 

2.4  Hailar  (IIIA) 

Station  HIA  is  the  northern  most  CDSN  station  located  on  the  Sino-Korean  Block  in  the 
Ne:  nenggu  Province  outside  Hailar  City  adjacent  to  the  Hailar  River.  The  seismometers  are 
located  in  a  vault  within  a  12  meter  long  tunnel  at  15  meters  depth  in  an  andesite  hill  (Peterson 
and  others — Kexin  and  others,  1987).  Topography  is  minimal  to  hilly  with  elevations  changes 
between  200-500  meters.  No  seismic  refraction  or  reflection  survey  results  in  the  vicinity  of  HIA 
were  found  in  the  literature.  Surface  wave  models  of  the  crust  and  upper  mantle  derived  from 
several  Kuril  earthquakes  recorded  by  Kimos  stations  WMQ,  KTA,  and  BTO  (Rui  et  al.  1981), 
provide  initial  estimates  beneath  HIA  shown  in  Figure  7.  Although  these  results  are  based  on 
one  azimuth  and  three  paths  Rui  et  al.  (1981)  infer  a  crustal  thickness  of  46  km  near  1 1 1 A  (model 

3)  with  an  average  velocity  of  6.4  km/s.  The  upper  mantle  P  velocity  is  between  7. 8-8.0  km/s 
and  all  models  indicate  a  step  in  velocity  at  approximately  23  km  depth.  Yuan  et  al.  (1986) 
shows  a  39  km  thick  crust  beneath  HIA  (Figure  2),  but  the  trend  of  increasing  Moho  depth  from 
station  MDJ  toward  HIA  is  consistent  in  both  studies. 

2.5  Mi  panjiam;  (MI). I) 

Station  MDJ  is  the  eastern  most  CDSN  station  situated  in  the  Sino-Korean  Block  of  the 
eastern  Heilongjiang  Province  approximately  300  km  northwest  of  the  Russian  port  Vladivostok. 
The  seismometers  are  installed  within  a  vault  toward  the  end  of  a  60  meter  long  tunnel  at  50 
meters  depth  in  granite  outside  the  city  of  Mujanjiang  (Peterson  and  others  Kexin  and  others, 
1987 ).  No  seismic  refraction  or  reflection  results  were  found  in  the  literature  for  the  MDJ  area. 
Surface  wave  results  given  by  Rui  et  al.  (1981)  shown  in  Figure  7  infer  a  42  km  thick  crust 
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(model  1)  with  a  slightly  lower  upper  mantle  velocity  ~7.9  km/s  beneath  station  MDJ.  As 
discussed  in  the  observations  for  station  HIA,  station  MDJ  also  shows  a  step  in  P  velocity  at  ~23 
km  depth.  Yuan  et  al.  (1986)  indicate  a  35  km  thick  crust  (Figure  2)  beneath  MDJ. 

2.6  URLIMQI  (YVMQ) 

Station  WMQ  is  located  between  the  Junggar  Basin  and  the  Tarim  Basin  along  the  eastern 
Tien  Shan  of  the  Xinjiang  Province.  The  sensors  are  installed  in  a  vault  at  6  meters  depth 
tunneled  into  a  sandstone  hill  on  the  outskirts  of  Urumqi  City  (Peterson  and  others —  Kexin  and 
others,  1987).  No  seismic  refraction  and/or  reflection  survey  results  in  the  vicinity  of  WMQ 
were  found  in  the  literature.  Yuan  et  al.  (1986)  infer  a  crustal  thickness  of  approximately  46-50 
km  in  the  vicinity  of  WMQ. 
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Siberian  Craton 


Figure  la.  Mercator  projection  of  China  showing  the  major  physiographic  features  and  CDSN  station  locations.  From  east 
to  west  these  stations  are:  (MDJ)  Mudanjiang,  (HIA)  Hailar,  (BJI)  Beijing,  (KMI)  Kunming,  (LZH)  Lanzhou,  and 
(WMQ)  Urumqi.  Figure  compiled  from  (Gilder  et  al.,  1991;  Hsu  et  al.,  1990;  and  Shangyou,  1991) 
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Figure  lb.  The  Provinces  of  China  (Hsu  et  al.  1990)  and  CDS  N  locations. 


Figure  2.  Contour  map  of  Moho  depth  in  kilometers  across  China,  from 
Yuan  et  al.  (1986). 
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Figure  4.  Refraction  profiles  with  respect  to  station  BJI  (top)  and  the  inferred 
velocity  model  (bottom)  given  by  Changquan  and  Shixu  (1986)  for  the  Zhenglanqi  (NW) 
to  Baigezhuang  (SE)  profile.  A  low  velocity  zone  in  the  mid-  to  lower  crust  is  also 
suggested  by  Shedlock  and  Roecker  (1987)  for  this  area. 
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Figure  5.  Tectonic  map  of  the  Yunnan-KMI-South  China  block  region  modified  from 
Baiji  and  Ning  (1987)  (top);  and  P-wave  velocity  model  from  Qizhong  et  al.  (1985) 
for  the  Simao-Malong  profile,  south  of  Malong  and  east  of  station  KM1. 
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Figure  6.  Depth  to  top  of  crust-mantle  boundary  around  Lanzhou  (LZH) 
(top)  from  Xingxin  (1988).  Distribution  of  Quaternary  active  faults 
and  Cenozoic  basins  (bottom)  from  Burchfiel  et  al.  (1991). 
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3.0  receiver  functions 

The  response  of  the  crust  and  upper-mantle  to  an  incoming  teleseismic  P  wave  may  be 
represented  with  a  convolutional  model: 

z(t)  =e,(t)  *  s(t)  *  i(t)  (1) 

where  z(t)  represents  the  vertical  ground  motion,  ez(t)  represents  the  vertical  component  response 
of  the  local  structure  to  an  incoming  plane  P  wave  and  includes  the  effects  of  laterally 
homogeneous  features  and  laterally  varying  structures,  s(t)  represents  the  apparent  source-time 
function  and  includes  both  near  source  and  lower  mantle  propagation  effects  and  i(t)  represents 
the  recording  instrument  response.  For  a  single  station  receiver  function  study  the  assumption 
of  a  laterally  homogeneous  structure  is  usually  necessary.  The  effects  of  lateral  heterogeneity 
and  dipping  structure  can  be  examined  through  the  variation  of  the  receiver  response  as  a 
function  of  the  incident  P  wave  azimuth  Langston  (1977a),  Owens  and  Crosson  (1988)  and  this 
study.  The  horizontal  radial  component  of  the  response  of  a  laterally  homogeneous  structure  to 
an  incoming  plane  P-wave  consists  of  the  direct  P  arrival  followed  by  P-to-S  converted  and 
reverberating  phases  and  zero  tangential  response. 

The  effects  of  the  source  and  lower  mantle  propagation  must  be  removed  to  isolate  the  local 
earth  structure.  Shown  in  Figure  8  is  one  of  several  characteristic  teleseismic  events  from  the 
Fiji  area  recorded  at  CDSN  station  BJI.  The  source  and  lower-mantle  propagation  effects  are 
common  to  all  three  components  of  motion,  and  with  matched  instruments,  a  receiver  function  is 
defined  as  the  time  series  resulting  from  the  deconvolution  of  the  vertical  component  of  motion 
from  the  horizontal  component  of  motion.  Below  we  describe  Langston’s  (1979)  source 
equalization  procedure  with  the  modification  of  (Ammon  1991).  Let  the  radial  component  of 
motion  be  represented  by  r(t),  the  vertical  by  z(t),  the  radial  receiver  function  by  h(t),  and  the 
corresponding  Fourier  transforms  by  R(co),  Z(to),  and  H(cn)  respectively.  The  time  domain 
receiver  function  is  defined  by 


r(l)  =  z(t)  *  h(t) 
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(2), 


such  that  in  the  frequency  domain  we  have 


Z(o» 


(3). 


To  reduce  the  effects  of  noise  and  numerical  instabilities  the  deconvolution  is  performed  in  the 
frequency  domain  using  a  minimum  allowable  denominator  described  by 


H(oi)  =  -(-0))-Z  ^)  G(o», 


(4) 


where 


4>(oi)  =  max  {Z(oi)  Z*(oi),  c  ■  max  (Z(oi)  Z’(0)) }  } 


(5). 


The  constant  c  represents  the  minimum  amplitude  allowed  in  the  denominator  of  (5),  *  indicates 
complex  conjugation  and 


G(oi)  (6) 

The  constant  £,  normalizes  the  gaussian  filter  to  unit  amplitude  in  the  time  domain  to  produce 
smooth  estimates  of  (4)  while  the  constant  “a”  determines  the  width  of  the  gaussian  filter.  A 
gaussian  of  2.5  passes  frequencies  up  to  ~1.2  Hz  while  a  gaussian  of  1.0  passes  frequencies  up  to 
~0.4  Hz.  Filling  spectral  holes  in  the  denominator  to  avoid  division  by  very  small  numbers 
changes  the  overall  amplitude  of  the  receiver  function.  Once  a  suitable  water-level  parameter  is 
chosen  the  true  receiver  function  amplitude  is  obtained  by  adding  the  following  constraint  to  the 
equalization.  The  deconvolution  of  the  vertical  component  from  itself  using  the  same  water-level 
parameter  as  in  the  horizontal  component  equalizations  produces  a  time  series  with  a  maximum 
amplitude  of  one  (Ammon,  1991).  True  amplitude  receiver  functions  preserve  the  relative 
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amplitude  of  the  converted  phases  to  the  direct  P  arrival.  Tne  receiver  function  ray  diagram  in 
Figure  9  shows  the  response  of  an  incident  plane  P-wave  after  propagation  across  layer  ‘h’ 
determined  with  the  technique  given  in  Kennet  (19X3).  The  use  of  a  plane  wave  approximation 
is  justified  because  a  teleseismic  P  wave  has  a  relatively  large  constant  phase  velocity.  For  a 
horizontally  layered  structure  each  P  to  S  conversion  has  the  same  horizontal  slowness.  The  true 
amplitude  synthetic  receiver  functions  computed  for  the  teleseismic  range  at  10°  intervals 
(Figure  9)  show  that  the  amplitude  falloff  :.nd  moveout  of  multiple  convened  phases  as 
epicentral  distance  increases  is  easily  observed.  These  fundamental  characteristics  are  observed 
in  the  CDSN  data  and  are  shown  in  Section  5. 

Of  primary  concern  is  the  contamination  of  a  radial  receiver  function  with  waves  generated 
by  the  interaction  of  the  P-wavefield  and  lateral  velocity  variations.  Differences  between  stacked 
azimuthal  clusters  may  be  attributed  in  part  to  variations  in  the  local  structure  as  a  function  of 
azimuth  and  in  part  by  contamination  with  off-azimuth  arrivals.  The  amount  of  scattering  is 
represented  by  the  observed  motion  on  the  tangential  receiver  function.  Unlike  interactive  trial- 
arid  error  modeling  the  inversion  scheme  is  unable  to  qualitatively  judge  the  importance  of  a 
given  arrival  in  a  receiver  function.  Manginoet  al.  (1992)  demonstrate  a  pitfall  in  interpretation 
of  a  receiver  function  inversion  in  the  presence  of  significant  scattering.  To  avoid  potentially 
biased  inversion  results  from  small  scale  (relative  to  wavelength)  heterogeneities  we  initially 
analyze  CDSN  data  using  a  gaussian  filter  of  1.0,  Longer  period  data  tend  to  smooth  over  the 
effects  of  small  scale  crustal  heterogeneities,  so  we  should  expect  to  obtain  simpler  and 
potentially  less  biased  solution  models.  When  multiple  stacking  suites  are  available  we  use  those 
with  the  best  signal-to-noi.se  ratio. 


August  10,  1990;  Fiji  Islands:  DELTA  =  85°,  BAZ  =120°  Depth  =  373  km 


Figure  8.  An  example  waveform  from  the  20  samples/second  broadband  data  stream 
of  seismograph  station  Baijatuan,  Beijing  (BJI)  one  of  6  CDSN  stations  providing  data 
to  the  Global  Digital  Seismic  Network.  Above  is  the  First  60  seconds  of  data  normalized 
to  the  vertical  component.  The  observed  ground  motion  corresponds  to  an  incident 
teleseismic  P-wave  from  the  Fiji  region. 


21 


Surface 


Direct  P 


Radial 

Receiver 

Function 


Figure  9.  Receiver  Function  Ray  Diagram  (top)  and  corresponding  receiver 
function  (middle).  Except  for  the  first  arrival  all  upper  case  letters  denote  down¬ 
going  waves  and  all  lower  case  letters  denote  up-going  waves  according 
to  Bath  and  Steffanson  (1966).  Shown  below  are  receiver  functions  from  the 
same  model  at  30°  to  90°  epicentral  distances  at  10°  intervals. 
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4.0  The  Data 

The  CDSN  became  operational  with  stations  BJI,  LZN  and  WMQ  during  mid-October 
19X6.  By  mid-March  1987  stations  KM1  and  HIA  were  ‘on-line’  and  toward  the  end  of  1987 
station  MDJ  was  working.  Most  stations  are  co-located  at  previously  established  fiducial  sites 
that  bear  the  same  name  (Peterson  and  others — Kexin  and  others,  1987).  Three  additional 
CDSN  stations  are  operational  (Qiongzhong  (QIZ),  Enshi  (ENH)  and  Sheshan  (SSE))  but  do  not 
routinely  report  data  to  the  international  seismological  community.  Deployment  of  new  stations 
are  planned  for  the  Lhasa,  Hotan  and  Golmud  regions  (Peterson  and  others — Kexin  and  others, 
1987).  All  CDSN  station  coordinates  are  listed  in  Table  7. 

In  this  report  we  examine  the  20  sps  BB/BH  STS-1  data  band.  The  CDSN  records  ground 
motion  at  sample  rates  of  [0.1,  1.0,  20,  40]  samples  per  second  (sps)  on  four  discrete  bands 
defined  respectively  as:  VLP  (very  long  period).  LP  (long  period),  BB,  BH  (broadband)  and  SP 
(short  period).  Three  types  of  seismometers  are  deployed.  The  BB,  LP  and  VLP  signal  bands  are 
recorded  with  a  triaxial  set  of  force-balance  Streckeisen  STS-1  seismometers.  The  SP  band  is 
recorded  with  China  Institute  of  Geophysics  DJ-l-SP  and  Geotech  KS-3600  borehole 
seismometers.  Station  clocks  are  synchronized  to  a  common  time  code  broadcast  by  station 
BPM  located  in  central  China.  Shown  in  Table  9  the  sensitivity  for  each  BB  data  channel  is 
nearly  identical  at  all  stations. 

All  ‘potentially  useful’  teleseismic  events  between  1986-1991  with  a  body  wave  magnitude 
of  5.8  or  greater  were  retrieved  from  the  archive.  A  non-potential  useful  event  is  defined  as  any 
teleseismic  event  recorded  at  a  given  station  within  the  coda  of  a  previous  local,  regional  or 
teleseismic  event.  Not  all  events  were  within  teleseismic  distance  of  all  stations  and  all 
potentially  useful  events  were  not  always  recorded  at  every  station.  Shown  in  Figure  10  are  the 
station  locations  with  the  number  of  teleseismic  events  (1986-1991)  recorded  by  each  station 
with  the  ‘potentially  useful’  constraint  included  in  the  tabulation.  In  general,  the  northeast  and 
southeast  quadrants  contain  the  most  data,  given  the  location  of  mainland  China  with  respect  to 
the  circum-Pacific  scismogenic  zones,  while  source  regions  for  the  southwest  and  northwest 
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quadrants  are  relatively  sparse.  Tables  1-6  contain  the  pertinent  source  parameters  as  listed  in 
the  Preliminary  Determination  of  Epicenter  Bulletins. 

Once  in-house  all  data  were  converted  into  SAC  format  (Seismic  Analysis  Code;  Tapley  et 
al.  1991)  and  are  shown  utilizing  SAC  graphics  subroutines.  For  reference  a  flow  chart 
summarizing  procedures  used  throughout  this  report  is  shown  Figure  11.  All  data  were  corrected 
for  sensitivity  and  sorted  into  azimuthal  quadrants.  Quadrants  that  contain  several  clusters  of 
events  were  binned  according  to  a  common  back  azimuth  and  epicentral  distance.  Throughout 
this  procedure  events  with  emergent  first  arrivals  are  discarded.  Individual  events  within  each 
bin  were  cut  into  60  second  record  lengths  with  30  seconds  of  pre-signal  noise.  Events  without 
30  seconds  of  pre-event  data  were  cosine  tapered  and  zero  padded.  Those  events  with  impulsive 
first  arrivals  are  selected  with  onset  picked  at  the  positive  peak  amplitude  first  arrival  on  the 
vertical  component.  The  radial  component  is  defined  as  the  positive  direction  from  the  source 
pointing  toward  the  receiver.  The  tangential  component  is  defined  as  90°  clockwise  from  the 
positive  radial  direction.  Each  event  was  source  equalized  with  the  procedure  discussed  in 
Section  3  and  those  events  within  each  bin  with  stable  deconvolutions  are  stacked.  The  variance 
of  the  stacked  data  provide  the  ±  1  standard  deviation  bounds  used  to  constrain  acceptable  model 
synthetics  in  the  following  sections  of  this  report. 
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Figuic  10.  Distribution  of  all  teleseismic  events  not  arriving  within  the  coda  of  a  local,  regional 
or  previous  teleseismic  event  with  respect  to  and  recorded  by  each  CDSN  station,  available  in 
the  archives  for  the  years  1986  until  1991. 


c 


Search  PDE’s; 
assemble  data  set. 


I 


) 


Access  archives  and  convert  various 
formats  into  SAC  formal. 


c 


Preliminary  data  processing 


J 


Process  all  traces  and  correct  for 
sensivitily.  Sort  into  bins  for  each 
station  for  each  azimuth. 


c 


Results  &  Discussion 


J 


Figure  1 1 .  Flowchart  of  processing  steps  used  in  this  report. 
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5.0  Receiver  Function  Observations 

As  shown  in  Figure  9  the  model  synthetics  clearly  show  that  the  overall  amplitude  falloff  and 
moveout  of  the  first  order  multiples  can  be  clearly  observed  over  tens  of  degrees.  The 
availability  of  several  stacks  within  a  narrow  azimuthal  range  provides  an  opportunity  to  observe 
these  characteristics  in  the  C1XSN  data.  In  Figures  12-  IK  and  the  following  sections  we  show  the 
response  and  discuss  our  observations  for  each  CDSN  station.  We  discuss  the  stacked  radial  and 
tangential  receiver  functions  with  the  best  signal-to-noise  ratio  for  quadrants  that  contain  several 
bins  and  we  also  compare  stacks  within  the  same  quadrant.  To  facilitate  comparison  of  arrivals 
common  to  all  a/imuths,  the  radial  and  tangential  receiver  functions  are  shown  normalized  to  the 
radial  component  within  each  quadrant  and  aligned  side  by  side.  Direct  comparison  between 
arrivals  on  the  radial  to  the  tangential  response  indicate  the  variation  in  time  of  scattered  energy 
to  the  energy  polarized  in  the  vertical-radial  plane.  As  a  rule  of  thumb  receiver  functions  sample 
a  lateral  distance  from  the  station  roughly  equal  to  the  horizontal  distance  traveled  by  the  deepest 
multiple  (typically  from  the  crust  mantle  boundary)  recorded  within  the  first  30s  after  the  direct  P 
arrival. 

5.1  BJI,  BKIJINt; 

The  distribution  of  data  and  the  5  data  bins  available  for  BJI  are  shown  in  Figures  A1-A6. 
Data  from  the  northeast  quadrant  consists  of  2  event  clusters  at  45°  and  58°  epicentral  distance 
with  an  azimuthal  separation  of  11°.  Southeast  of  BJI  bins  at  74°and  85°  each  contain  several 
events.  Three  receiver  functions  were  sufficient  to  form  a  single  bin  for  data  arriving  from  the 
west. 

Shown  in  Figure  12  are  the  BJI  receiver  functions.  Common  to  all  azimuths  is  a  large 
coherent  pulse  of  energy  on  the  radial  response  arriving  at  5-5. 3s  after  the  direct  arrival.  After  7- 
Xs  both  radial  and  tangential  responses  are  significantly  different.  The  BJI-NE  radial  receiver 
function  appears  relatively  simple  with  a  large  pulse  of  energy  at  5s  after  the  direct  arrival 
followed  by  lower  energy  arrivals.  Tangential  motion  from  ()-3s  show  a  moderate  level  of 


energy  with  a  positive  first  motion  followed  by  energy  close  to  but  not  exceeding  the  amplitude 
on  radial.  The  BJI-SH  response  is  complicated.  From  7s  to  about  18s  the  radial  response 
contains  several  large  amplitude  arrivals  while  the  tangential  response  has  a  small  negative 
motion  at  onset  followed  by  arrivals  that  exceed  the  radial  after  6s.  The  pronounced  scattered 
energy  on  the  tangential  component  prohibits  a  confident  interpretation  of  this  azimuth  in  terms 
of  a  one-dimensional  structure.  The  BJ1-NW  radial  response  has  three  coherent  pulses  of  energy. 
The  pulse  at  5.3s  is  followed  by  a  large  positive  motion  at  15s  and  an  equally  prominent  negative 
motion  at  20s.  The  tangential  response  shows  moderate  energy  from  onset  to  10s  but  does  not 
exceed  the  radial.  After  10s  the  tangential  response  is  minimal  relative  to  the  radial. 

One  interpretation  of  the  arrival  at  5-5. 3s  common  to  all  azimuths  is  a  Ps  conversion  at  the 
crust-mantle  boundary.  Furthermore,  the  similarity  of  this  arrival  about  the  station  suggests  that 
the  region  of  the  Moho  structure  sampled  by  this  Ps  phase  is  laterally  consistent  about  the  station. 
Motion  on  the  NW  response  between  14-20s  is  consistent  with  PpPhS  and  PpShS  multiples  from 
the  Moho.  From  these  observations  only  the  northeast  and  northwest  azimuths  shows  promise  for 
obtaining  reliable  one-dimensional  model  estimates  of  the  crustal  structure. 

5.2  K.Ml,  KiiNMiNt; 

Station  KMI  has  8  bins  shown  in  Figures  A7-A14.  Data  from  the  northeast  contains  2  event 
clusters  at  approximately  63°  and  75°.  The  southeast  quadrant  consists  of  3  bins  at  approximately 
56°,  75c,  and  88°.  From  source  regions  to  the  west  two  bins  are  possible  with  the  northwest 
quadrant  containing  several  events. 

Shown  in  Figure  13  tire  the  KMI  receiver  functions.  In  general  arrivals  from  the  east  are 
similar  and  arrivals  from  the  west  are  similar,  but  between  east  and  west  there  are  significant 
differences  in  amplitude  and  number  of  arrivals.  The  northeast  tangential  response  has  large 
energy  immediately  from  onset  until  about  20s.  The  radial  response  contains  large  energy  pulses 
at  5-5. 2s  and  between  l()s  to  20s  that  is  in  phase  with  motion  on  the  tangential.  Large  amplitude 
tangential  arrivals  preclude  interpretation  of  this  azimuth  in  terms  of  a  one-dimensional  response. 


The  KMI-SE  response  is  also  complex.  Tangential  motion  through  the  initial  15s  is  significant 
but  low  in  amplitude  compared  to  the  radial  component.  Between  5-8s  the  southeast  radial 
response  contains  the  largest  energy  arrival  observed  at  this  station,  followed  by  a  second  packet 
of  (multiple?)  energy  between  1  l-l 8s.  The  KM1  western  receiver  functions  are  lower  amplitude 
after  about  1  Os  compared  to  the  receiver  functions  from  the  east  and  motion  on  the  radial 
components  appear  to  be  polarized  in  the  vertical-radial  plane.  These  observations  are  consistent 
with  a  one-dimensional  response  but  does  not  preclude  an  alternate  structure. 

5.3  LZH,  Lanzhou 

Station  LZH  contains  7  event  clusters  shown  in  Figures  A15-A22.  The  northeast  quadrant 
contains  2  bins  at  55°  and  68°  with  an  azimuthal  separation  of  14°.  The  southeast  quadrant 
contains  4  bins  at  46°,  62°,  77°  and  93°  consisting  of  5,  4,  7  and  4  events  respectively.  This 
coverage  is  the  best  of  any  station  quadrant.  West  of  LZH  a  single  azimuthal  cluster  consisting 
of  4  events  is  possible. 

Shown  in  Figure  14  are  the  LZH  receiver  functions.  Although  the  radial  and  tangential 
response  are  different  for  each  backazimuth,  all  radial  components  show  a  moderately  coherent 
pulse  of  energy  between  6-7s  after  the  direct  arrival.  The  NF  tangential  response  contains  large 
energy  from  onset  to  about  12.  The  SE  radial  response  contains  a  large  energy  packet  between  2- 
8s  after  the  direct  arrival.  Between  8- 14s  the  tangential  motion  is  large  but  does  not  exceed  the 
radial  response.  Shown  in  Figure  15  spanning  almost  50°  great  circle  arc  are  all  LZH-SF  radial 
and  tangential  receiver  functions.  Differences  in  each  bin  suggest  a  difference  in  crustal  structure 
over  the  35°  azimuthal  range.  As  the  ray  parameter  decreases  the  radial  response  approaches  the 
energy  level  on  the  tangential  response.  The  amplitude  fall  off  of  the  5-7s  pulse  and  moveout  of 
the  pulse  between  13- 14s  (inset)  are  consistent  with  a  conversion  from  a  common  interface 
beneath  LZH. 
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The  LZH-W  response  is  relatively  simple  compared  to  the  east.  The  packet  of  energy  at  5-8s 
on  the  radial  component  is  the  only  clearly  distinguishable  arrival.  The  western  tangential 
component  has  the  lowest  energy  of  all  azimuths  and  does  not  exceed  the  radial  response. 

5.4  HIA,  HAILAR 

Station  HIA  contains  6  azimuthal  clusters  shown  in  Figures  A23-A30.  The  northeast 
quadrant  consists  of  one  bin  at  50°  while  the  SE  quadrant  is  well  sampled  by  3  bins  at  55°,  62° 
and  73°.  Data  arriving  from  the  west  can  be  grouped  into  2  bins  at  34°  and  50°. 

Shown  in  Figure  16  are  the  HIA  radial  and  tangential  receiver  functions.  In  general,  from  the 
direct  P  arrival  to  about  8-9s  the  overall  response  is  very  similar.  The  pulse  between  4-5s 
followed  by  a  smaller  arrival  between  6-7s  is  common  to  all  rauial  components.  Only  the  SE 
radial  response  contains  an  arrival  between  13-1%  that  is  consistent  with  a  multiple  from  the 
same  interface  generating  the  4-5s  arrival.  Tangential  energy  is  minimal  with  respect  to  the 
radial  for  all  azimuths.  AH  HIA-SE  stacks  are  shown  in  Figure  17.  The  strong  similarity  of  the 
radial  response  for  each  bin  indicates  the  crustal  structure  is  similar  over  the  azimuthal  range 
sampled.  The  decreasing  overall  amplitude  and  the  moveout  of  a  suspected  PpPhS  multiple 
(Figure  17)  as  a  function  of  epicentral  distance  is  consistent  with  model  synthetics  shown  in 
Figure  9.  Moveout  of  a  multiple  from  a  common  interface  can  provide  additional  constraint  to 
the  velocity  structure  beneath  HIA. 

5.5  MDJ,  Mi  danjian<; 

Station  MDJ  was  the  last  CDSN  station  designated  to  provide  data  to  the  GDSN  to  become 
fully  operational  in  1987.  Subsequently  MDJ  data  have  been  sporadic  in  the  archives.  Shown  in 
F  igures  A30-A33  are  individual  MDJ  receiver  functions.  In  particular  events  9212  (SE)  and 
9213  (SW)  are  consistent  with  a  one-dimensional  response.  Overall,  with  only  a  handful  of 
events  distributed  about  the  station  we  believe  that  an  analysis  for  the  receiver  structure  beneath 
MDJ  is  unwarranted  at  this  time. 


5.6  YVMQ,  IJRUMQI 

The  VVMQ  data  contains  5  possible  bins  shown  in  Figure  A35-A40.  The  northeast  quadrant 
consists  of  a  single  well  populated  bin  at  63°.  The  southeast  quadrant  contains  two  bins  at  63° 
and  75°,  separated  by  21°.  Data  arriving  from  western  sources  yield  one  bin  at  50°  in  the 
southwest  and  a  second  bin  at  3K°  in  the  northwest  quadrant. 

Shown  in  Figure  18  are  the  WMQ  receiver  functions.  The  overall  difference  in  response 
indicates  a  laterally  varying  velocity  structure  beneath  the  station.  Motion  on  the  NE  tangential 
is  the  largest  while  motion  on  the  SE  tangential  is  the  least  for  all  azimuths.  Close  examination 
of  the  radial  response  reveals  an  offset  of  the  direct  arrival  for  all  azimuths.  This  observation 
coupled  with  a  significant  variation  in  polarity  within  the  initial  (Ms  of  the  tangential  response 
suggest  the  presence  of  a  dir>p  ig  structure  beneath  WMQ.  Owens  and  Crosson  (1988)  show  that 
the  time  lag  of  the  '\r  i  arrival  on  the  radial  response  is  related  to  the  traveltime  of  a  Ps 
conversion  propagating  through  the  up-dip  or  down-dip  direction,  provided  that  velocity 
increases  with  depth.  A  Ps  conversion  traveling  up-dip  has  the  largest  delay  because  of  a  longer 
travel  path  through  a  slower  medium.  Langston  (1977)  shows  that  the  initial  motion  on  the 
tangential  response  is  reversed  across  a  line  parallel  to  the  dip  direction. 

In  light  of  these  previous  studies  to  constrain  the  suspected  dipping  structure  we  measured 
the  time  lag  for  all  radial  receiver  functions  and  examine  all  tangential  receiver  functions  for 
polarity  reversals  as  a  function  of  backazimuth.  Shown  in  Figure  19a  are  the  observed  time-lags 
measured  from  all  stable  WMQ  radial  receiver  functions.  A  consistent  pattern  with  an  average 
delay  of  ().8s  is  clearly  present  between  310°  to  15°.  Figure  19b  shows  the  tangential  response 
about  WMQ.  Arrows  indicate  the  backazimuth  range  over  which  we  can  constrain  a  line  parallel 
to  the  dip  direction.  The  time  delays  of  the  direct  arrival  plus  the  polarity  reversals  are  consistent 
with  a  north-northwest  dipping  structure  beneath  the  station.  The  effects  of  this  feature  are  most 
likely  diverting  converted  shear  wave  energy  out  of  the  vertical-radial  plane  and  therefore  could 
lead  to  biased  inversion  results. 
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Figure  13.  KMI  Receiver  Functions 
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Figure  14.  LZH  Receiver  Functions 


LZH-SE  Response  with  increasing  epicentrol  distance 
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Figure  15 


Figure  16.  HIA  Receiver  Functions 
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Figure  18.  WMQ  Receiver  Functions 
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Time  delay  (seconds) 


Figure  19a.  Shown  above  are  the  time  lags  for  all  first  arrivals  on 
the  radial  component.  A  clear  time  delay  as  a  function  of  azimuth 
suggests  a  dipping  structure  toward  the  NW  of  the  station. 


J9 


WMQ  tangential  receiver  functions 


Figure  19b.  Distribution  of  tangential  receiver  functions  by  backazimuth  at  WMQ. 
Polarity  reversal  of  the  First  arrival  occurs  along  a  line  parallel  to  the  dip  direction 
(Langston  1977).  Arrows  indicate  the  backazimuthal  range  where  the  polarity  of 
the  tangential  reverses. 
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6.0  Receiver  Function  Inversion 

The  inversion  procedure  used  in  this  report  and  the  notation  summarized  below  are  consistent 
with  those  presented  Ammon  et  al.  (1990)  who  give  a  rigorous  description  of  the  method.  We 
highlight  here  those  points  most  important  for  understanding  and  interpreting  our  results.  The 
non-linear  relationship  between  the  velocity  structure  and  the  waveform  is  approximated  by  a 
Taylor  series  expansion  about  an  initial  model,  m0: 

dj  =  Fjfnt]  =  Fjfm,,]  +  (D,8m)j  +  OlSm^l  j=  1.2,3,  ..N  (7) 

where  dj  represents  the  observed  data,  m  is  the  true  model,  parentheses  denote  the  inner  product, 
and  D  is  a  matrix  containing  the  partial  derivatives  of  the  waveform  with  respect  to  changes  in 
the  layer  velocities  in  ni„.  The  operator  F  represents  the  nonlinear  relationship  between  the 
model  and  the  receiver  function,  and  it  includes  the  calculation  of  a  plane-wave  synthetic 
seismogram  and  the  source  equalization  procedure  described  in  Section  3.0.  Neglecting  the 
higher  order  terms  (7)  becomes 

(D,5m)j  =  Fj[m]  -  Fj[m„]  (S) 

The  right  hand  side  of  (X)  is  the  waveform  residual  vector  and  (X)  may  be  solved  for  8m  using 
standard  least-squares  techniques  described  by  Wiggins  (1972)  and  Aki  and  Richards  (19X0). 
Adding  (D,  m„)  to  both  sides  of  (X)  enables  solution  directly  for  m  to  produce: 

(D,m)j  =  Fj[m]  -  F^m,,]  +  (D,m„)j 

Equation  (9)  is  a  matrix  equation  with  the  only  unknown  being  m.  By  solving  directly  for  the 
model  m,  a  smoothness  constraint  can  be  applied  during  the  solution  of  (9).  We  refer  the  reader 
to  Ammon  et  al.  (1990)  on  the  topic  of  applying  smoothness  constraints  to  the  solution  of  (9). 
The  matrix  D  is  calculated  by  employing  Randall's  (19X9)  efficient  differential  seismogram 


41 


calculation  technique.  Compressional  wave  velocities  are  determined  using  a  constant  Poisson's 
ratio  of  0.25  and  the  layer  density  is  adjusted  by  using  the  equation  p  =  0.32  Vp  +  0.77 
(Berteussen  1977). 

In  practice,  once  an  initial  model  is  obtained  from  previous  work,  and  from  forward  modeling, 
we  run  a  set  of  inversions  with  a  range  of  smoothness  values  in  order  to  minimize  both  model 
roughness  and  RMS  error  between  the  model  synthetic  and  the  stacked  data.  With  the  pseudo¬ 
monte-carlo  technique  the  initial  model  is  randomly  perturbed  into  24  different  starting  models 
by  adding  a  random  cubic  polynomial  function  as  well  as  a  random  component  for  each  layer. 
The  non-uniqueness  of  the  problem  is  examined  by  inverting  the  receiver  function  beginning 
with  many  different  starting  models  (Ammon  et  al.  1990).  The  first  inversion  suite  has  a  random 
cubic  perturbation  of  0.5  kin/s  and  the  second  has  inversion  suite  has  a  1.0  km/s  perturbation, 
which  allows  a  larger  range  of  starting  models  to  test  the  sensitivity  of  solution  models  on  the 
initial  model.  Shown  in  Figure  20  are  the  starting  models  for  station  BJI.  A  total  of  48  starting 
models  are  inverted  and  all  model  synthetics  that  do  not  fit  the  most  coherent  arrivals  within  the 
variance  of  the  stacked  data  are  discarded.  The  range  of  remaining  solution  models  reflects  the 
non-uniqueness  of  the  particular  problem. 


h'l 


Vp  in  km/sec 


\  6  8 


5  20.  Range  of  starting  models  for  station  BJI. 
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7.0  PRELIMINARY  RESULTS 

In  this  section  we  present  preliminary  inversion  results  for  each  CDSN  station.  Our  goal  here 
is  to  determine  model  estimates  of  the  receiver  structure  by  forcing  only  the  most  prominent 
phases  in  the  synthetic  receiver  functions  to  be  within  the  ±1  standard  deviation  bounds  obtained 
from  the  variance  of  the  stacked  data.  By  relaxing  the  smoothness  constraint  as  well  as  shifting 
the  layer  thickness  and  upper  mantle  velocity  we  may  better  fit  the  smaller  amplitude  arrivals  in 
the  stacked  data.  However  in  some  cases  the  lower  amplitude  arrivals  on  the  radial  are  close  to 
the  energy  level  present  on  the  tangential.  Therefore  caution  must  be  used  when  matching  such 
energy  solely  with  P  to  S  conversions  generated  within  the  vertical-radial  plane.  For  those 
stations  where  the  tangential  response  is  small  the  initial  model  solutions  should  provide  a  good 
estimate  of  the  velocity  structure  beneath  the  station. 

7.1  BJI 

Shown  in  Figure  21  are  preliminary  model  results  for  station  BJI.  The  starting  model  was 
estimated  from  previous  refraction  models  near  the  station.  In  general  the  crust  west  of  station 
BJI  appears  relatively  simple  in  structure.  The  upper  6  km  is  characterized  by  a  strong  positive 
gradient.  Between  10  to  33  km  depth  the  average  crustal  P  wave  velocity  is  6  km/s,  and  a  small 
low  velocity  zone  is  present  at  about  15  km  depth.  Between  33-39  km  depth  a  strong  positive 
gradient  is  present,  and  the  top  of  the  upper  mantle  is  at  39  km  depth. 

7.2  K  M  I 

Shown  in  Figure  23  are  preliminary  model  estimates  for  KMI.  No  models  have  been 
selected  based  on  synthetics  that  fit  the  data  within  one  standard  deviation.  The  starting  model 
was  estimated  from  previous  refraction  models  near  the  station.  The  wide  range  of  crustal 
velocities  suggests  significant  differences  exist  between  the  starting  model  and  the  structure 
necessary  to  fit  the  data.  Shifting  the  layer  thickness,  upper  mantle  velocity,  and  decreasing  the 
smoothness  constraint  may  narrow  the  range.  The  variability  in  solution  models  at  KMI  is  not 


surprising  given  the  strong  waveform  variability  of  the  radial  receiver  function  and  the  large 
tangential  energy  as  a  function  of  azimuth. 

7.3  L7J1 

Shown  in  Figure  24  are  preliminary  model  estimates  for  LZH.  The  synthetics  are  strongly 
dependent  on  changes  in  the  near-surface  velocity.  Model  synthetics  from  a  high  velocity 
shallow  crust  fall  on  and  above  the  upper  resolution  bound,  while  synthetics  from  lower  velocity 
shallow  crustal  models  fall  on  and  below  the  lower  resolution  bound.  Shifting  layer  thickness, 
upper  mantle  velocity,  and  decreasing  the  smoothness  constraint  may  narrow  the  range  of 
possible  crustal  models. 

7.4  HIA 

Preliminary  results  for  station  HIA  are  shown  in  Figure  25.  The  model  synthetics  provide  a 
reasonably  good  fit  to  the  data.  The  starting  model  was  estimated  from  surface  wave  estimates 
near  HIA.  The  upper-crust  has  a  shallow  jump  in  velocity  followed  by  a  negative  gradient  to  -7.5 
km  depth.  The  mid-crust  has  a  relatively  constant  velocity  of  -6.0  km/sec  to  a  depth  of  28  km, 
and  overlies  a  transitional  crust-mantle  boundary  between  30-37  km  depth.  The  top  of  the  upper 
mantle  is  at  37-38  km  depth. 

7.5  WMQ 

Shown  in  Figure  26  are  preliminary  results  at  station  WMQ  for  data  arriving  in  the  up-dip 
direction.  The  upper  crust  shows  a  strong  positive  gradient  to  about  6-8  km  depth,  which  is 
necessary  to  fit  the  large  time  lag  measured  for  this  azimuth.  The  absence  of  a  prominent 
converted  phases  results  in  a  smoothly  varying  velocity  structure  with  increasing  depth. 
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8.0  DISCUSSION 

The  main  advantage  of  the  receiver  function  technique  is  that  it  can  be  used  to  obtain 
estimates  of  the  velocity  structure  beneath  a  seismic  station  without  resorting  to  specialized  field 
studies  around  the  station.  Furthermore,  the  method  is  able  to  resolve  the  velocity  structure  of 
the  crust  and  upper  mantle  within  only  a  few  tens  of  kilometers  beneath  a  seismic  station  as 
compared  to  seismic  refraction  techniques  which  typically  estimate  the  crust  and  upper  mantle 
structure  averaged  over  50-100  km.  However,  a  tradeoff  between  receiver  function  and 
refraction  methods  arises  because  refraction  experiments  can  record  waveforms  with  frequencies 
up  to  25-30  Hz  while  this  receiver  function  study  examines  frequencies  up  to  about  0.4  Hz. 
Thus,  while  low-pass  receiver  functions  tend  to  smooth  over  the  fine  details  of  the  crustal 
structure,  they  are  able  to  constrain  depth-to-interface  and  velocity  values  for  the  major  velocity 
layers  beneath  a  station. 

8.1  Implications  for  the  Monitoring  of  Nuclear  Test  Ban  Treaties 

Determination  of  the  velocity  structure  beneath  remote  seismic  stations  has  important 
applications  for  the  monitoring  of  existing  nuclear  test  ban  treaties  and  clandestine  testing. 
Accurate  event  locations  require  a  detailed  knowledge  of  the  seismic  velocity  structure  of  the 
earth  (Blandford,  1982),  something  which  is  not  well  known  in  potential  areas  where  clandestine 
testing  could  take  place.  Also,  accurate  knowledge  of  Earth  structure  beneath  the  recording 
station  can  be  an  important  factor  in  refining  yield  estimates  of  nuclear  explosions  (Bache,  1982, 
Hansen  et  al.,  1990).  Some  of  the  variability  in  amplitude  and  waveshape  measurements  of  the 
body  waves  (typically  the  P  waves)  from  underground  nuclear  explosions  is  due  to  scattering  and 
body-wave  conversions  beneath  the  receiving  stations  (e.g.,  Bannister,  1990).  The  receiver 
function  technique  is  useful  in  addressing  the  effect  of  this  particular  problem  at  a  given  seismic 
station. 

Shown  in  Figure  26  are  great  circle  paths  from  the  September  14,  1988  JVE2  nuclear 
explosion  at  the  Soviet  test  site  in  Semipalatinsk  (SEM)  to  the  five  CDSN  stations  that  recorded 
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this  event.  The  source-receiver  distances  are:  WMQ  -  950  km;  LZH  -  2530  km;  HIA  -  2926  km; 
BJI  -  3106  km;  and  KMI  -  3432  km.  The  effects  of  the  velocity  structure  beneath  these  CDSN 
stations  on  the  waveforms  from  this  event  are  explored  in  Figures  27a, b.  Figure  27a  shows 
Filtered  observations  (vertical  and  radial)  for  the  P  wave  and  some  of  the  P  coda.  Synthetic 
waveforms  for  these  stations  are  shown  in  Figure  27b.  The  synthetics  were  computed  by  taking 
the  Mueller-Murphy  source  pulse  shape  for  an  118  kt  explosion  in  granite  and  convolving  it 
with  the  appropriate  instrument  response  and  the  vertical  and  radial  responses  computed  from 
one  of  the  mode's  determined  in  this  study  for  each  station.  The  synthetics  are  quite  simple. 
They  do  not  include  important  upper  mantle  triplications  which  must  be  present  in  most  of  the 
seismograms,  nor  do  they  include  free  surface  reflections  (pP),  spall,  source  structure  or  tectonic 
release.  The  purpose  of  the  synthetics  is  to  document  how  much  variability  in  the  P  waveforms 
can  be  expected  due  to  the  differences  in  the  receiver  structures  found  beneath  the  CDSN 
stations.  Figure  27b  shows  the  same  major  effects  documented  by  Bannister  et  al.  (1990)  in  the 
P  waves  from  NORESS.  Motions  on  the  radial  component  are  much  more  affected  by  the 
estimated  receiver  structures  than  the  vertical  component.  However,  even  for  the  vertical 
component,  the  width  of  the  second  peak  for  the  P  wave  (the  so-called  "e"  peak)  varies  from 
station  to  station.  Furthermore,  the  relative  amplitude  of  the  "c"  peak  to  the  first  (or  "a")  peak 
differs  somewhat  among  the  stations  for  the  vertical  P  component.  The  radial  ground  motions 
show  variations  among  the  width  of  the  "c"  peaks  as  well  as  in  the  waveforms  about  5  seconds 
after  the  P-wave.  The  latter,  hints  of  which  may  be  observed  in  the  data,  reflect  differences  in 
the  crustal  structure  from  station  to  station.  Thus,  as  concluded  by  Bannister  et  al.  (1990)  P-to-S 
scattering  does  have  a  significant  effect  upon  the  P  coda  whereas  P-to-P  scattering  does  not. 
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e  26.  Mercator  projection  of  great  circle  paths  between  Semipalatinsk  and  CDSN  stations  (HIA)  Hailar,  (BJI)  Beijing, 
(KMI)  Kunming,  (LZH)  Lanzhou,  and  (WMQ)  Urumqi. 
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Figure  27a.  JVE2  vertical  (left)  and  radial  (data)  as  recorded  at 
all  CDSN  stations  bandpass  filtered  between  0.2 — 1.1  Hz. 
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Figure  27b.  Synthetic  vertical  and  radial  seismograms  obtained  from  a  characteristic 
preliminary  model  from  each  station,  convolved  with  a  JVE2  source  time  function  and 
STS-1  instrument  response. 
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Table  1  Station  BJI 


Event 

Latitude0 

Longitude0 

Mb 

BAZ° 

Delta0 

Location 

NE 

1 

0121 

58.840  N 

156.858  W 

6.1 

39 

55 

Alaska  ** 

2 

1052 

58.419  N 

1 75.445  W 

6.3 

42 

46 

Bering  Sea 

3 

7005 

52.448  N 

169.381  W 

6.1 

49 

50 

Fox  Island  * 

4 

7049 

51.298  N 

179.279  W 

6.2 

53 

45 

Andreanof  Islands  * 

5 

7321 

58.586  N 

143.270  W 

6.6 

35 

61 

Gulf  of  Alaska  ** 

6 

8038 

50.785  N 

173.465  E 

6.2 

55 

40 

Aleutian  Islands  * 

7 

8066 

56.953  N 

143.032  W 

6.8 

37 

62 

Gulf  of  Alaska  ** 

8 

8C85 

62.154  N 

124.182  W 

6.1 

26 

67 

Canada 

9 

9247 

55.543  N 

156.835  W 

6.5 

42 

56 

Alaska  ** 

10 

9280 

51.314  N 

1 79.028  W 

6.1 

53 

45 

Aleutian  Islands  * 

SE 

11 

0177 

22.015  S 

179.473  W 

6.0 

123 

86 

Fiji* 

12 

0222 

19.085  S 

177.385  W 

6.0 

120 

85 

Fiji* 

13 

1008 

18.016  S 

173.660  W 

6.1 

116 

87 

Tonga 

14 

6303 

21.702  S 

176.616  W 

6.4 

120 

87 

Fiji* 

15 

7041 

19.489  S 

177.456  W 

6.2 

119 

85 

Fiji  * 

16 

7078 

20.397  S 

176.134  W 

5.9 

119 

87 

Fiji  * 

17 

7119 

19.013  S 

1 77.736  W 

5.9 

119 

85 

Fiji* 

18 

9070 

17.766  S 

174.761  W 

6.4 

119 

85 

Tonga* 

19 

0208 

15.355  S 

167.464  E 

6.4 

128 

72 

Vanuatu  “ 

20 

0224 

19.435  S 

169.132  E 

6.3 

129 

77 

Vanuatu  ** 

21 

0271 

13.559  S 

167.079  E 

6.0 

127 

71 

Vanuatu  ** 

22 

7042 

15.834  S 

167.355  E 

5.9 

128 

73 

Vanuatu 

23 

8157 

15.397  S 

167.578  E 

6.0 

127 

72 

Vanuatu  ** 

24 

9096 

19.306  S 

169.002  E 

6.1 

129 

76 

Vanuatu 

25 

9268 

20.355  S 

169.277  E 

6.1 

129 

77 

Vanuatu  ** 

NW 

26 

0151 

45.811  N 

26.769  E 

6.1 

307 

62 

Romania  * 

27 

0310 

28.251  N 

55.462  E 

6.2 

276 

50 

Iran* 

28 

1119 

42.489  N 

43/647  E 

6.2 

297 

53 

Western  Caucasus 

29 

7125 

36.480  N 

70.673  E 

5.8 

279 

35 

Hindu  Kush  Region 

30 

7352 

28.191  N 

56.677  E 

5.8 

275 

49 

Iran 

31 

8342 

40.987  N 

44.485  E 

6.2 

296 

53 

Turkey-USSR  border 

Event  location  modifications  correspond  to  the  following  azimuthal  bins:  *BJI-45-NE, 
**BJI-58-NE;  *BJI-86-SE,  **BJI-74-SE.  See  Appendix  1  for  individual  events. 
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Table  2  Station  KMI 


Event 

Latitude0 

Longitude' 

3  Mb 

BAZ° 

Delta0 

Location 

NE 

1 

7049 

51.298  N 

179.279  W 

6.2 

43 

63 

Andreanof  1st. 

2 

7058 

53.470  N 

167.291  W 

6.2 

39 

70 

Fox  Island  ** 

3 

7080 

52.056  N 

177.547  W 

6.0 

42 

64 

Andreanof  1st.  * 

4 

7081 

51.594  N 

173.574  W 

5.9 

42 

66 

Andreanof  1st.  * 

5 

7172 

54.211  N 

162.601  W 

6.2 

37 

72 

Alaska  Peninsula  " 

6 

8038 

50.785  N 

173.465  E 

6.2 

44 

58 

Aleutian  Islands  * 

7 

8320 

52.109  N 

171.103  W 

5.9 

41 

68 

Fox  Island 

8 

8066 

56.953  N 

143.032  W 

6.8 

30 

81 

Gulf  of  Alaska 

9 

8085 

62.154  N 

124.182  W 

6.1 

20 

85 

Canada 

10 

9139 

54.305  N 

165.574  W 

6.1 

38 

70 

Aleutian  Islands 

11 

9167 

31.807  N 

137.982  E 

5.9 

32 

75 

Honshu 

12 

9247 

55.543  N 

156.835  W 

6.5 

35 

75 

Alaska  " 

13 

9280 

51.314  N 

179.028  W 

6.1 

43 

63 

Aleutian  1st.  * 

14 

0121 

58.840  N 

156.858  W 

6.1 

32 

74 

Alaska  ** 

SE 

15 

7041 

19.489  S 

1 77.456  W 

6.2 

111 

89 

Fiji A 

16 

7042 

15.834  S 

167.355  E 

5.9 

115 

74 

Vanuatu  '* 

17 

7285 

7.288  S 

154.371  E 

6.3 

115 

59 

Solomon  1st. 

18 

7338 

5.777  S 

154.555  E 

5.9 

114 

58 

Solomon  1st.  * 

19 

8022 

19.847  S 

133.803  E 

6.2 

143 

53 

Northern  Australia  * 

20 

8070 

20.917  S 

178.645  W 

6.1 

113 

89 

Fiji 

21 

8137 

13.941  S 

166.335  E 

6.0 

114 

73 

Vanuatu 

22 

8157 

15.397  S 

167.578  E 

6.0 

115 

74 

Vanuatu  ** 

23 

8184 

14.278  S 

167.180  E 

5.9 

114 

73 

Vanuatu  ** 

24 

8209 

13.1 12  S 

167.051  E 

5.9 

113 

73 

Vanuatu  ** 

25 

8205 

6.526  S 

152.779  E 

6.7 

116 

57 

New  Britain  Region 

26 

8205a 

22.127  S 

174.900  E 

5.9 

117 

84 

Loyalty  1st A 

27 

8219 

7.136  S 

151.057  E 

5.9 

117 

56 

New  Britain  Region 

28 

9035 

4.625  S 

153.066  E 

6.1 

113 

57 

New  Ireland  Region  * 

29 

9070 

17.766  S 

174.761  W 

6.4 

113 

89 

Tonga A 

30 

9268 

6.526  S 

152.779  E 

6.7 

118 

78 

Vanuatu  ** 

31 

0177 

22.015  S 

1 79.473  W 

6.0 

114 

88 

Fiji A 

32 

0208 

22.015  S 

1 79.473  W 

6.0 

115 

74 

Vanuatu  ** 

33 

0224 

19.435  S 

169.132  E 

6.3 

117 

78 

Vanuatu  “ 

34 

0271 

13.559  S 

167.097  E 

6.0 

113 

73 

Vanuatu  “ 

35 

0364 

5.097  S 

150.967  E 

6.6 

115 

55 

New  Britain  Region  * 

SW 

36 

9069 

13.702  S 

34.420  E 

6.2 

248 

76 

Malawi  * 

37 

9232 

1 1.766  N 

41.942  E 

5.8 

268 

58 

Ethiopia  * 

38 

9233 

1 1.874  N 

41 .870  E 

6.3 

268 

58 

Ethiopia  * 

39 

0362 

14.875  S 

66.777  E 

6.0 

225 

53 

Mid- Indian  Rise 

NW 

40 

8342 

40.987  N 

44.485  E 

6.2 

303 

50 

Turkey-USSR  border  * 

41 

8219 

36.461  N 

71.043  E 

6.1 

300 

29 

Afghn.-USSR  border 

42 

9259 

40.337  N 

51.534  E 

6.4 

302 

45 

Caspian  Sea  * 

43 

9260 

42.489  N 

43.647  E 

6.2 

302 

44 

Caspian  Sea  * 

44 

0150 

45.341  N 

26.688  E 

6.7 

310 

62 

Romania 

45 

0310 

28.251  N 

55.462  E 

6.2 

285 

42 

Iran  * 

46 

1119 

42.489  N 

43/647  E 

6.2 

305 

51 

Western  Caucasus  * 

Event  location  modifications  correspond  to  the  following  azimuthal  bins:  *KMI-63-NE 

**KMI-73-NE; 

*KMI-56-SE,  **KMI-75-SE,  AKMI-88-SE.  See  Appendix  1  for 

individual  events. 


Table  3  Station  LZH 


Event 

Latitude0 

Lonqitude0 

Mb 

BAZ° 

Delta0 

Location 

NE 

1 

7049 

51.298  N 

1 79.279  W 

6.2 

48 

55 

Andreanof  Isl.* 

2 

7080 

52.056  N 

1 77.547  W 

6.0 

46 

55 

Andreanof  Isl.* 

3 

7081 

51.594  N 

173.574  W 

5.9 

46 

58 

Andreanof  Isl. 

4 

7172 

54.211  N 

162.601  W 

6.2 

40 

63 

Alaska  Peninsula 

5 

7279 

52.956  N 

159.972  E 

6.1 

48 

42 

Kamchatka 

6 

7321 

58.586  N 

143.270  W 

6.6 

30 

70 

Gulf  of  Alaska  ** 

7 

8038 

50.785  N 

173.465  E 

6.2 

50 

50 

Aleutian  Islands 

8 

8047 

51.564  N 

175.041  E 

5.9 

46 

57 

Rat  Isl.  * 

9 

8066 

56.953  N 

143.032  W 

6.8 

32 

71 

Gulf  of  Alaska  ** 

10 

8067 

41.669  N 

152.222  E 

5.9 

66 

37 

Northern  Pacific 

11 

8085 

62.154  N 

124.182  W 

6.1 

21 

74 

Canada 

12 

8320 

52.109  N 

171.103  W 

5.9 

45 

59 

Fox  Island 

13 

9022 

41.806  N 

144.282  E 

6.0 

61 

31 

Hokkaido 

14 

9065 

35.545  N 

140.444  E 

5.9 

80 

29 

Honshu 

15 

9101 

49.488  N 

159.185  E 

6.3 

53 

41 

Kuril  Islands 

16 

9139 

54.305  N 

165.574  W 

6.1 

41 

61 

Aleutian  Islands 

17 

9144 

56. 1 77  N 

164.264  E 

5.9 

43 

44 

Komandorsky  Islands 

18 

9167 

31.807  N 

137.982  E 

5.9 

34 

66 

Honsu 

19 

9247 

55.543  N 

156.835  W 

6.5 

37 

65 

Alaska  ** 

20 

9280 

51.314  N 

1 79.028  W 

6.1 

48 

55 

Aleutians  * 

21 

9282 

51.780  N 

171.869  E 

6.0 

49 

49 

Aleutians 

22 

9355 

45.364  N 

150.108  E 

5.9 

60 

35 

Kuril  Isl 

23 

0121 

58.840  N 

156.858  W 

6.1 

34 

64 

Alaska  ** 

SE 

24 

6303 

21.702  S 

176.616  W 

6.4 

113 

94 

Fiji  * 

25 

7003 

14.998  S 

167.929  E 

6.0 

92 

60 

Vanuatu 

26 

7041 

19.489  S 

1 77.456  W 

6.2 

112 

92 

Fiji A 

27 

7078 

20.397  S 

1 76.134  W 

5.9 

112 

94 

Fiji A 

28 

7168 

5.577  S 

130.791  E 

6.6 

143 

48 

Banda  Sea  * 

29 

8022 

19.847  S 

133.803  E 

6.2 

147 

62 

Northern  Australia  ** 

30 

8116 

7.791  S 

158.255  E 

6.1 

119 

67 

Solomon  Islands  ** 

31 

8124 

22.744  S 

170.278  E 

5.9 

121 

85 

Loyalty  Isl.  Region 

32 

8137 

13.941  S 

166.335  E 

6.0 

117 

77 

Vanuatu 

33 

8151 

7.501  S 

128.325  E 

6.5 

147 

49 

Banda  Sea  * 

34 

8157 

15.397  S 

167.578  E 

6.0 

118 

78 

Vanuatu  - 

35 

8184 

14.278  S 

167.180  E 

5.9 

117 

77 

Vanuatu- 

36 

8205 

6.526  S 

152.779  E 

6.7 

122 

62 

New  Britain  Region 

37 

8209 

13.112  S 

167.051  E 

5.9 

116 

77 

Vanuatu  - 

38 

9010 

3.162  S 

130.556  E 

5.9 

141 

46 

Ceram  * 

40 

9067 

1.031  N 

126.189  E 

5.9 

146 

41 

Molucca  Passage 

41 

9195 

8.081  S 

125.129  E 

6.4 

151 

48 

Timor' 

42 

9257 

1.644  N 

127.322  E 

6.0 

142 

40 

Haimahera ' 

43 

9268 

6.526  S 

152.779  E 

6.7 

120 

83 

Vanuatu  - 

44 

9291 

10.155  S 

161 .065  E 

6.1 

118 

70 

Solomon  Islands 

45 

9300 

11.022  S 

162.350  E 

6.1 

118 

72 

Solomon  Islands  - 

46 

0177 

22.015  S 

179.473  W 

6.0 

115 

92 

Fiji  A 

47 

0208 

22.015  S 

179.473  W 

6.0 

118 

78 

Vanuatu  - 

48 

0271 

13.559  S 

167.097  E 

6.0 

117 

72 

Vanuatu  - 

49 

0364 

5.097  S 

150.967  E 

6.0 

122 

60 

New  Britain  Region  ** 

50 

1003 

7  .148  S 

148.490  E 

6.0 

126 

60 

East  Papua  New  Guinea  ** 

Table  3  Station  LZH  (continued) 


Event 

Latitude0 

Lonqitude0 

Mb 

BAZ° 

Delta0 

Location 

51 

7352 

28.191  N 

56.677  E 

5.8 

272 

40 

Iran 

52 

8094 

4.687  N 

94.419  E 

5.9 

197 

32 

Northern  Summatera 

53 

8342 

40.987  N 

44.485  E 

6.2 

295 

46 

Turkey-USSR  border 

54 

9259 

40.337  N 

51.534  E 

6.4 

292 

40 

Caspian  Sea  * 

55 

9260 

42.489  N 

43/647  E 

6.2 

292 

40 

Caspian  Sea 

56 

0310 

28.251  N 

55.462  E 

6.2 

273 

41 

Iran* 

57 

0362 

14.875  S 

66.777  E 

6.0 

221 

61 

Mid  Indian  Rise 

53 

0150 

45.841  N 

26.688  E 

6.7 

305 

56 

Romania 

5S 

0151 

45.811  N 

26.769  E 

6.1 

305 

56 

Romaina 

60 

0310 

28.251  N 

55.462  E 

6.2 

273 

41 

Iran  * 

61 

1119 

42.489  N 

43/647  E 

6.2 

297 

46 

Western  Caucasus  * 

Event  location  modifications  correspond  to  the  following  azimuthal  bins:  *LZH-57-NE, 
**LZH-68-NE;  *LZH-46-SE,  **LZH-62-SE,  -LZH-77-SE,  ALZH-93-SE.  See  Appendix 
1  for  individual  events. 


Table  4  Station  HIA 


Event 

Latitude0 

Lonqitude0 

Mb 

BAZ° 

Delta0 

Location 

NE 

1 

8066 

56.953  N 

143.032  W 

6.8 

42 

54 

Gulf  of  Alaska  * 

2 

8085 

62.154  N 

124.182  W 

6.1 

29 

57 

Canada  * 

3 

9280 

51.314  N 

179.028  W 

6.1 

62 

38 

Aleutian  Islands  * 

4 

8320 

52.109  N 

171.103  W 

5.9 

58 

42 

Fox  Island 

SE 

5 

9195 

8.018  S 

125.129  E 

6.4 

173 

57 

Timor ' 

6 

9212 

8.048  S 

121.384  E 

6.3 

178 

57 

Flores  Isl  Region  * 

7 

8151 

7.501  S 

128.325  E 

6.5 

169 

57 

Banda  Sea  * 

8 

9010 

3.162  S 

130.556  E 

5.9 

166 

53 

Ceram  * 

9 

7168 

5.577  S 

130.791  E 

6.6 

166 

55 

Banda  Sea  * 

10 

7178 

2.164  S 

138.170  E 

5.7 

156 

53 

West  Irian  * 

11 

9213 

4.511  S 

139.022  E 

6.0 

156 

56 

West  Irian 

12 

0208 

15.355  S 

167.464  E 

6.4 

132 

76 

Vanuatu  ** 

13 

0271 

13.559  S 

167.079  E 

6.0 

132 

75 

Vanuatu  ** 

14 

8137 

13.941  S 

166.335  E 

6.0 

133 

75 

Vanuatu  ** 

15 

8157 

15.397  S 

167.578  E 

6.0 

132 

77 

Vanuatu  ** 

16 

8184 

14.278  S 

167.180  E 

5.9 

132 

75 

Vanuatu  ** 

17 

9035 

10.447  S 

161.372  E 

6.0 

135 

69 

Solomon  Isl. 

18 

9291 

10.155  S 

161 .063  E 

6.1 

136 

69 

Solomon  Isl. 

19 

8116 

7.791  S 

158.255  E 

6.1 

137 

66 

Solomon  Isl. "  A 

20 

0364 

5.097  S 

150.967  E 

6.6 

143 

60 

New  Britain  Region A 

21 

1003 

7.148  S 

148.490  E 

6.0 

147 

61 

East  Papua,  New  Guinea 

22 

7285 

7.288  S 

154.371  E 

6.3 

141 

63 

Solomon  Isl. A 

23 

8190 

6.290  S 

154.667  E 

5.8 

140 

63 

Solomon  Isl. 

24 

8205 

6.526  S 

152.779  E 

6.7 

142 

62 

New  Britain  Region  A 

25 

8219 

7.136  S 

151.057  E 

5.9 

144 

62 

New  Britain  Region A 

26 

1008 

18.016  S 

173.660  W 

6.1 

119 

89 

Tonga 

27 

9096 

19.306  S 

169.002  E 

6.1 

133 

81 

Vanuatu 

sw 

28 

0310 

28.251  N 

55.462  E 

6.2 

271 

52 

Iran 

29 

1005 

23.478  N 

95.983  E 

6.3 

224 

32 

Burma* 

30 

7125 

36.480  N 

70.673  E 

5.8 

268 

37 

Hindu  Kush  Region  * 

31 

7352 

28.191  N 

56.677  E 

5.8 

270 

52 

Iran 

32 

8233 

26.755  N 

86.616  E 

6.4 

240 

34 

Nepal -India  * 

33 

9232 

11.766  N 

41.942  E 

5.8 

268 

73 

Ethiopia 

34 

0362 

14.875  S 

66.777  E 

6.0 

231 

79 

Mid-Indian  Rise 

NW 

35 

1119 

42.489  N 

43/647  E 

6.2 

293 

51 

Western  Caucasus  * 

36 

8342 

40.987  N 

44.485  E 

6.2 

291 

51 

Turkey-USSR  border  * 

37 

9260 

40.203  N 

51.749  E 

6.1 

286 

47 

Caspian  Sea  * 

38 

0150 

45.841  N 

26.688  E 

6.7 

305 

58 

Romania 

39 

0151 

45.811  N 

26.769  E 

6.1 

305 

58 

Romania 

Event  location  modifications  correspond  to  the  following  azimuthal  bins:  *HIA-55-SE, 
AHIA-62-SE,  **HIA-73-SE.  See  Appendix  1  for  individual  events. 


Table  5  Station  MDJ 


Event 

Latitude0 

Lonaitude0 

Mb 

BAZ° 

Delta0 

Location 

NE 

1 

9139 

54.305  N 

165.574  W 

6.1 

52 

42 

Fox  Islands 

2 

9247 

55.543  N 

156.835  W 

6.5 

48 

46 

Alaska 

SE 

3 

9134 

30.523  S 

1 78.414  W 

5.9 

137 

88 

Kermadec 

4 

9135 

9.803  S 

159.531  E 

5.9 

145 

61 

Solomon  Isl 

5 

9140 

30.508  S 

178.270  W 

5.7 

137 

88 

Kermadec 

6 

9213 

4.511  S 

139.022  E 

6.0 

167 

50 

West  Irian 

7 

9268 

20.355  S 

169.277  E 

6.1 

141 

74 

Vanuatu 

8 

9291 

10.155  S 

161 .063  E 

6.1 

144 

62 

Solomon  Isl 

9 

8022 

19.847  S 

133.803  E 

6.2 

175 

64 

Northern  Australia 

SW 

10 

9163 

21.861  N 

89.763  E 

6.1 

248 

40 

Bangladesh 

11 

9212 

8.048  S 

121.348  E 

6.3 

191 

53 

Flores  Isl.  Region 

12 

9257 

1.644  N 

127.322  E 

6.0 

183 

43 

Halmahara 

NW 

13 

9233 

1 1.874  N 

41.870  E 

6.3 

276 

80 

Ethiopia 

14 

9259 

40.337  N 

51.543  E 

6.4 

295 

55 

Caspian  Sea 

15 

9260 

40.203  N 

51.749  E 

6.1 

294 

55 

Caspian  Sea 

See  Appendix  1  for  individual  receiver  functions. 

Table  6  Station  WMQ 


Event 

Latitude0 

Lonaitude0 

Mb 

BAZ° 

Delta0 

Location 

NE 

1 

7005 

52.448  N 

169.381  W 

6.1 

41 

63 

Fox  Island  * 

2 

7058 

53.470  N 

167.291  W 

6.2 

40 

63 

Fox  Island  * 

3 

7080 

52.056  N 

1 77.547  W 

5.9 

45 

59 

Andreanof  Isl. 

4 

7081 

51.594  N 

173.574  W 

5.9 

44 

61 

Andreanof  Isl.  * 

5 

7172 

54.211  N 

162.601  W 

6.2 

37 

65 

Alaska  Peninsula  * 

6 

8038 

50.785  N 

173.465  E 

6.2 

50 

55 

Aleutian  Islands 

7 

8066 

56.953  N 

143.032  W 

6.8 

26 

70 

Gulf  of  Alaska 

8 

8085 

62.154  N 

124.182  W 

6.1 

15 

71 

Canada 

9 

9065 

35.545  N 

140.444  E 

5.9 

83 

40 

Honshu 

10 

9082 

33.788  N 

141.395  E 

5.6 

84 

42 

Honshu 

11 

9167 

31.807  N 

137.982  E 

5.9 

30 

66 

Honshu 

12 

9280 

51.314  N 

179.028  W 

6.1 

46 

59 

Aleutian  Isl  * 

13 

0121 

58.840  N 

156.858  W 

6.1 

31 

64 

Alaska* 

14 

9247 

55.543  N 

156.835  W 

6.6 

33 

66 

Alaska* 

SE 

15 

7168 

5.577  S 

130.791  E 

6.6 

130 

63 

Banda  Sea  * 

16 

7178 

2.164  S 

138.170  E 

5.7 

121 

64 

West  Irian  * 

17 

7338 

5.777  S 

154.555  E 

5.9 

110 

77 

Solomon  Islands  ** 

18 

8094 

4.687  N 

94.419  E 

5.9 

169 

39 

Northern  Sumatera 

19 

8116 

7.791  S 

158.255  E 

6.1 

109 

81 

Solomon  Islands 

20 

8151 

7.501  S 

128.325  E 

6.5 

133 

63 

Banda  Sea  * 

21 

8157 

15.397  S 

167.578  E 

6.0 

107 

93 

Vanuatu 

22 

8184 

14.278  S 

167.108  E 

5.9 

107 

92 

Vanuatu 

23 

8205 

6.526  S 

152.779  E 

6.7 

112 

76 

New  Britain  Region  ** 

24 

8207 

6.081  S 

133.667  E 

6.5 

127 

65 

Aroe  isl  * 

25 

8209 

13.112  S 

167.051  E 

5.9 

106 

91 

Vanuatu 

26 

9010 

3.162  S 

130.556  E 

5.9 

128 

61 

Ceram  * 

27 

9059 

2.288  N 

129.916  E 

5.7 

127 

54 

Molucca  Passage 

28 

9067 

1.031  N 

126.189  E 

5.9 

132 

54 

Molucca  Passage 

29 

9067a 

1.031  N 

126.189  E 

5.9 

130 

55 

Molucca  Passage 

30 

9121 

4.200  S 

101.366  E 

5.5 

161 

49 

Southern  Sumatera 

31 

9135 

9.803  S 

159.531  E 

5.9 

110 

84 

Solomon  Islands 

32 

9144 

3.748  S 

159.564  E 

5.6 

105 

79 

Sulawesi 

33 

9165 

12.875  N 

143.351  E 

5.5 

105 

56 

South  of  Mariana  Isl 

34 

9195 

8.081  S 

125.129  E 

6.4 

136 

62 

Timor* 

35 

9212 

8.048  S 

121.384  E 

6.3 

140 

60 

Flores  Isl  Region  * 

36 

9213 

4.511  S 

139.022  E 

6.0 

122 

67 

West  Irian* 

37 

9291 

10.155  S 

161.063  E 

6.1 

108 

85 

Solomon  Isl. 

38 

9300 

1 1.022  S 

162.305  E 

6.1 

108 

86 

Solomon  Isl. 

39 

0364 

5.097  S 

150.967  E 

6.6 

112 

74 

New  Britain  Region  ** 

40 

1003 

7.148  S 

148.490  E 

6.0 

116 

74 

East  Papua  New  Guinea  ** 

SW 

41 

9068 

13.710  S 

34.381  E 

5.8 

234 

75 

Malawi 

42 

9069 

13.702  S 

34.420  E 

6.2 

234 

75 

Malawi 

43 

9232 

11.766  N 

41.942  E 

5.8 

245 

50 

Ethiopia  * 

44 

9233 

1 1.874  N 

41.870  E 

6.3 

245 

50 

Ethiopia  * 

45 

0063 

28.925  S 

66.331  E 

6.8 

199 

75 

Pakistan 

46 

0362 

14.875  S 

66.777  E 

6.0 

203 

61 

Mid  Indian  Rise 

65 


Table  6  Station  WMQ  (continued) 


Event 

Latitude0 

Lonqitude0 

Mb 

BAZ° 

Delta0 

Location 

NW  47 

8342 

40.987  N 

44.185  E 

6.2 

280 

32 

Turkey-USSR  border  * 

48 

9160 

71.432  N 

4.371  W 

5.6 

335 

49 

Jan  Mayen  Isl. 

49 

9177 

39.112  N 

28.242  W 

5.7 

314 

79 

Azores 

50 

9359 

60.008  N 

73.445  W 

6.2 

350 

75 

Northern  Quebec 

51 

0150 

45.841  N 

26.769  E 

6.7 

294 

42 

Romania  * 

52 

0151 

45.81 1  N 

26.769  E 

6.1 

294 

42 

Romania  * 

Event  location  modifications  correspond  to  azimuthal  bins:  *WMQ-63-SE,  **WMQ-75- 
SE.  See  Appendix  1  for  individual  events. 
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Table  7  CDSN  STATION  LOCATIONS 


BJl  tat:  40.040277°  N 

Lon:  116.175000°  E 
Elev:  +43.0  meters 
Baijatuan,  Beijing 

HIA  Lat:  49.266666° N 

Lon:  119.7416666°  E 
Elev:  +610.0  meters 
Hailar,  Neimenggu  Province 

KMI  Lat:  25.148055°  N 

Lon:  102.747222°  E 
Elev:  +1952.0  meters 
Kunming,  Yunnan  Province 

LZH  Lat:  36.086666°  N 

Lon:  103.84444°  E 
Elev:  +1560.0  meters 
Lanzhou,  Gansu  Province 

MDJ  Lat:  44.616388°  N 

Lon:  129.591944°  E 
Elev:  +250.0  meters 
Mudanjiang,  Heilongjiang  Province 

WMQ  Lat:  43.821 111°N 

Lon:  87.695°  E 
Elev:  +970.0  meters 
Urumqi,  Xinjiang  Province 

QIZ  Lat:  19.029444°  N 
Lon:  109.843333°  E 
Elev:  +230  meters 
Qiongzhong,  Hainan  Island 

SSE  Lat:  31.095555°  N 

Lon:  121.186666°  E 
Elev:  +10  meters 
Sheshan 


ENH  Lat:  46.266666°  N 
Lon:  109.497500°  E 
Elev:  +487  meters 
Enshi,  Hubie  Province 


Table  8.  Sensitivity  of  each  CDSN  BB  Channel 


BJI.BBZ 

BJI.BBN 

BJI.BBE 

7.200E+09 

6.890E+09 

6.910E+09 

HIA.BBZ 

HIA.BBN 

HIA.BBN 

7.120E+09 

7.020E+09 

7.020E+09 

KM1.BBZ 

KMI.BBN 

KMI.BBE 

6.810E+09 

6.960E+09 

7.020E+09 

LZH.BBZ 

LZH.BBN 

LZH.BBE 

7.470E+09 

6.830E+09 

6.900E+09 

MDJ.BBZ 

MDJ.BBN 

MDJ.BBE 

7.340E+09 

7.020E+09 

7.020E+09 

WMQ.BBZ 

WMQ.BBN 

WMQ.BBE 

7.160E409 

6.870E+09 

6.850E+09 

(counts/meter/second) 
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APPENDIX 


In  the  following  style  of  file  identification  each  radial  receiver  function  corresponds  to  each 
single  or  double  asterisk,  hat,  or  tilde  appended  to  a  given  event  in  Tables  1-5.  Event  indexes  are 
defined  by  calender  year  and  Julian  day.  For  example  0241  .bji.eqr  corresponds  to  1990-Julian 
day  241-station  BJl-radial  receiver  function.  Shown  below  each  bin  is  the  stack  and  the  bin 
filename  includes  the  average  epicentral  distance. 


Figure  A- 1 
Figure  A-2 
Figure  A-3 
Figure  A-4 
Figure  A-5 
Figure  A-6 


BJ1  event  distribution. 

BJI  Northeast  45°  data  and  stack. 
BJI  Northeast  58°  data  and  stack. 
BJI  Southeast  74°  data  and  stack. 
BJI  Southeast  85°  data  and  stack. 
BJI  Northwest  56°data  and  stack 


Figure  A-7 
Figure  A-8 
Figure  A-9 
Figure  A- 10 
Figure  A- 1 1 
Figure  A- 12 
Figure  A- 13 
Figure  A- 14 


KMI  event  distribution. 

KM1  Northeast  63°  data  and  stack. 
KMI  Northeast  75°  data  and  stack. 
KMI  Southeast  56°  data  and  stack. 
KMI  Southeast  75°  data  and  stack. 
KMI  Southeast  88°  data  and  stack. 
KMI  Southwest  64°data  and  stack. 
KMI  Northwest  46°  data  and  stack. 


Figure  A- 15 
Figure  A-16 
Figure  A- 17 
Figure  A- 18 
Figure  A-19 
Figure  A-20 
Figure  A-21 
Figure  A-22 


LZH  event  distribution. 

LZH  Northeast  55°  data  and  stack. 
LZH  Northeast  68°  data  and  stack. 
LZH  Southeast  46°  data  and  stack. 
LZH  Southeast  62°  data  and  stack. 
LZH  Southeast  77°  data  and  stack. 
LZH  Southeast  93°  data  and  stack. 
LZH  West  43°  data  and  stack. 


Figuie  A-23 
Figure  A-24 
Figure  A-25 
Figure  A-26 
Figure  A-27 
Figure  A-28 
Figure  A-29 


HIA  event  distribution. 

HIA  Northeast  50°  data  and  stack. 
HIA  Southeast  55°  data  and  stack. 
HIA  Southeast  62°  data  and  stack. 
HIA  Southeast  73°  data  and  stack. 
HIA  Southwest  34°  data  and  stack. 
HIA  Northwest  50°  data  and  stack. 


Figure  A-30 
Figure  A-31 
Figure  A-32 
Figure  A-33 
Figure  A-34 


MDJ  data  distribution 
MDJ  Northeast  data. 
MDJ  Southeast  data. 
MDJ  Southwest  data. 
MDJ  Northwest  data. 


Figure  A-35  WMQ  Event  Distribution. 

Figure  A-36  WMQ  Northeast  63°  data  and  stack. 

Figure  A-37  WMQ  Southeast  75°  data  and  stack. 

Figure  A-38  WMQ  Southeast  63°  data  and  stack. 

Figure  A-39  WMQ  Southwest  50°  data  and  stack. 

Figure  A-40  WMQ  Northwest  38°  data  and  stack. 
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BJI  Event  Distribution 


S 


Figure  A-l.  Shaded  ellipses  correspond  to  azimuthal  bins  shown  for 
each  quadrant. 


70 


BJI-45-NE 


BJI-58-NE 
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Figure  A-4 


BJI-56-NW 


Figure  A-6 


73 


KMI  Event  Distribution 


F  igure  A-7.  Shaded  ellipses  correspond  to  azimuthal  bins  shown  for 
each  quadrant. 
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Figure  A-8 
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Figure  A- 10 
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Figure  A- 1  I 
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KU1-88-SE 
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Figure  A- 14 


78 


L2H-55-NE 


LZH-6B-NE 


is 


•JO  -20  -10  0  10  20 

Seconds 


Figure  A- 18 
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Figure  A- 19 


81 


LZH-77-SE 


LZH-93-SC 
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Figure  A-22 
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Figure  A-27 
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HIA-34-SW 


Figure  A-28 


HA-50-N* 
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Figure  A-31 
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Figure  A-32 


89 


UQJ  southwestern  quodront 


Figure  A-33 
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Figure  A-34 
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WMQ  Event  Distribution 
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Figure  A-35.  Shaded  ellipses  correspond  to  azimuthal  bins  shown 
for  each  quadrant. 
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Figure  A-39 
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WUQ-38-NW 


Figure  A -40 
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